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(— I Abstract 

Cd In this paper we consider the spectral and nonlinear stability of periodic traveling 

H wave solutions of a generalized Kuramoto-Sivashinsky equation. In particular, we re- 

'— ^ solve the long-standing question of nonlinear modulational stability by demonstrating 

^^ that spectrally stable waves are nonlinearly stable when subject to small localized (in- 

<^ tegrable) perturbations. Our analysis is based upon detailed estimates of the linearized 

l/~j solution operator, which are complicated by the fact that the (necessarily essential) 

0^ spectrum of the associated linearization intersects the imaginary axis at the origin. We 

r carry out a numerical Evans function study of the spectral problem and find bands of 

; spectrally stable periodic traveling waves, in close agreement with previous numerical 

rn studies of Frisch-She-Thual, Bar-Nepomnyashchy, Chang-Demekhin-Kopelevich, and 

others carried out by other techniques. We also compare predictions of the associated 

Whitham modulation equations, which formally describe the dynamics of weak large 

scale perturbations of a periodic wave train, with numerical time evolution studies, 

demonstrating their effectiveness at a practical level. For the reader's convenience, we 
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include in an appendix the corresponding treatment of the Swift-Hohenberg equation, 
a nonconservative counterpart of the generaHzed Kuramoto-Sivashinsky equation for 
which the nonhnear stabihty analysis is considerably simpler, together with numerical 
Evans function analyses extending spectral stability analyses of Mielke and Schneider. 
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1 Introduction 

Localized coherent structures such as solitary waves play an essential role as elementary 
processes in nonlinear phenomena. Examples of this are multi-bump solutions in reaction 
diffusion equations, which are constructed by piecing together well-separated solitary waves 
[E], or the limiting case of infinite, periodic wave trains. A similar situation occurs in non- 
linear dispersive media described by a KdV equation where exact multi-bump and periodic 
solutions exist. In this paper, we consider periodic solutions of an unstable dissipative- 
dispersive nonlinear equation, namely a generalized Kuramoto-Sivashinsky (gKS) equation 

(1.1) ut + -fd^u + ed^u + 5dlu + d^f{u) = 0, 7, <^ > 0, 

where f{u) is an appropriate nonlinearity and e, 7 G M are arbitrary constants with 7 > 0. 

2 

In the case f{u) = "^^ equation (1.1) is a canonical model for pattern formation that has 
been used to describe, variously, plasma instabilities, flame front propagation, turbulence in 
reaction-diffusion systems and nonlinear waves in fluid mechanics [SI, S2, K, KT, CD, PSU]. 

Equation (1.1) may be derived formally either from shallow water equations [YY] or 
from the full Navier-Stokes system [W] for < 7 = 5 ^ 1. Here 5 measures the deviation 
of the Reynolds number from the critical Reynolds number above which large scale weak 
perturbations are spectrally unstable. For this latter application, what we have in mind is 
the description of chaotic motions in thin film flows down an incline [CD] . Indeed, periodic 
traveling waves are some of the few simple solutions in the attractor for the classic (e = 0) 
Kuramoto-Sivashinsky equation, a generic equation for chaotic dynamics, and there is now 
a huge literature on these solutions (and their bifurcations, in particular period doubling 
cascades) and their stability; see [EST, KE, K, CD]. As e increases, the set of stable periodic 
waves, and presumably also their basin of attraction appears (numerically) to enlarge [CDK, 
BaN], until, in the |(7, (^)| — ?■ limit, they and other approximate superpositions of solitary 
waves appear to dominate asymptotic behavior [CD, PSU, BJRZ, BJNRZ2]. 

Since (5 > it is easily seen via Fourier analysis that all constant solutions of (1.1) 
are unstable, from which it follows that all asymptotically constant solutions (such as the 
solitary waves) are also unstable. Nevertheless, one can still construct multi-bump solutions 
to (gKS) on asymptotically large time 0{5~^) by gluing together solitary waves, provided 
that the distance between them is not too large [PSU]. One possible interpretation of this 
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is that there exist stable periodic wave trains nearby the sohtary wave. Indeed, it has been 
known, almost since the introduction of the classical Kuramoto-Sivashinsky equation (1.1) 
(e = 0) in 1975, that there exists a spectrally stable band of periodic solutions in parameter 
space; see for example the numerical studies in [CKTR, FST]. These stable periodic wave 
trains may be heuristically viewed as a superposition of infinitely many well separated 
solitary waves. In [EMR], the existence of such a band of stable periodic traveling waves 
was justified for the equation (1.1) with periodic boundary conditions and in the singular 
KdV limit 1(7,(5)1 ^0. 

Although numerical time-evolution experiments suggest that these spectrally stable 
waves are nonlinearly stable as well (see [CD]), up to now this conjecture had not been 
rigorously verified. In this paper, relying heavily on the recent infusion of new tools in 
[JZl, JZN, BJRZ, BJNRZl, BJNRZ2] in the context of general conservation laws and 
the St. Venant (shallow water) equations, we prove the result, previously announced in 
[BJNRZ3], that spectral modulational stability of periodic solutions of (1.1), defined in the 
standard sense, implies linear and nonlinear modulational stability to small localized (inte- 
grable) perturbations; that is, a localized perturbation of a periodic traveling wave converges 
to a periodic traveling wave that is modulated in phase. The first such nonlinear result for 
any version of (1.1), this closes in particular the 35-year old open question of nonlinear 
stability of spectrally stable periodic waves of the classical Kuramoto-Sivashinsky equation 
(e = 0) found numerically in [FST]. 

With these improvements in nonlinear theory, we find this also an opportune moment 
to make a definitive discussion of the generalized Kuramoto-Sivashnisky equation (and 
Swift-Hohenberg equations) in terms of existence, nonlinear theory, and numerical spectral 
stability studies, all three, across all parameters, both connecting to and greatly generalizing 
the variety of prior works [FST, KE, K, CD, Ml, Sc, CDK, BaN]. We thus carry out also 
a numerical analysis of the spectrum in order to check the spectral assumptions made in 
our main theorem. Since the spectrum always contains the origin, due to both translational 
invariance and the conservative form of the equations, this stability is at best marginal, a 
circumstance that greatly complicates both numerical and analytical stability analyses. 

Our numerical approach is based on complementary tools; namely Hills method and the 
Evans function. On the one hand, we use SpectrUW numerical software [CDKK] based 
on Hills method, which is a Galerkin method of approximation, in order to obtain a good 
overview of location of the spectrum: the periodic coefficients and eigenvectors are expanded 
using Fourier series, and then a frequency cutoff is used to reduce the problem to finding 
eigenvalues of a finite dimensional matrix. It is known that Hills method converges faster 
than any algebraic order [JZ5]; moreover, in practice, it gives quickly a reliable global 
qualitative picture. However, the associated error bounds are of abstract nature, with 
coefficients whose size is not a priori guaranteed. Further, near the critical zone around 
A = 0, the resolution of this method is not in practice sufficient to guess at stability, let 
alone obtain satisfactory numerical verification. 

Thus, in order to get more reliable pictures near the origin and guarantee the spectral 
stability of periodic wave, we use on the other hand an Evans function approach, computing 
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a winding number to prove that there is or is not unstable spectrum in the part of the 
unstable (positive real part) complex half plane excluding a small neighborhood of the origin, 
then using Cauchy's integral formula to determine the Taylor expansions of the spectral 
curves passing through the origin. This method, though cumbersome for approximating 
global spectrum, is excellent for excluding the existence of spectra on a given region, and 
comes with error bounds that can in a practical sense be prescribed via the tolerance of the 
Runga-Kutta 4-5 scheme used to evaluate the Evans function by solution of an appropriate 
ODE; see [Br, Zl, BJNRZl] for general discussion of convergence of Evans function methods. 
Furthermore, under generic assumptions, the numerical protocol introduced in Section 2.1.3 
below detects sideband stability and instability of the underlying periodic wave train without 
the need of lengthy spectral perturbation expansion calculations, thus adding what we 
believe is a valuable new method to the numerical toolbox for analyzing the spectral stability 
of periodic wave trains. For relations between Hill's method, Fredholm determinants, and 
the periodic Evans function of Gardner, see [JZ5, Z2]. 

In order to validate our numerical method, we check benchmarks from the Kuramoto- 
Sivashinsky equation {e = 0, f{u) = v? jl) and the Swift Hohenberg equation. We obtain 
very good agreement with existing numerical works. In particular, we recover and extend 
stability boundaries found numerically for KS in [EST] and analytically for Swift-Hohenberg 
in [Ml, Sc]. Then we perform an all-parameters study linking in particular the Kuramoto- 
Sivashinsky equation e— )'07 = 5=lto KdV equation £ = 0,7 = 5— )'0 by homotopy, 
again obtaining excellent agreement with existing results of [CDK, BaN]. A rigorous proof 
of spectral stability of periodic traveling wave solutions of (1.1), in the spirit of [Ml, Sc], or 
a numerical proof (interval arithmetic) would be a very interesting direction for future work, 
particularly in an interesting parameter regime such as the canonical KDV limit | (7, (J)| — > 
studied in [PSU] . In this regard, we point to the recent singular perturbation analysis of this 
limit in [BaN] and [JNRZ4], the second relying on and completing the first, reducing the 
problem of rigorous validation of stability/instability as |(7, 5)| — > to computer- verified 
evaluation of the signs of certain elliptic integrals associated with KdV. 

As we will see, our main theorem indicates that the long time dynamics of a localized 
perturbation of a given periodic traveling wave of (1.1) is governed by a space-time depen- 
dent phase modulation ^(x,t), which in turn satisfies in some sense a set of modulation 
equations. Moreover, it is well known that, in the neighborhood of the origin, the spectrum 
of the linearization about a given periodic wave train is strongly related to some local well- 
posedness properties of the associated linearized Whitham system; this set of modulation 
equations may be derived via a nonlinear optics expansion (WKB), and formally govern the 
evolution of weak large scale perturbations of the underlying wave train. 

In [EST], the authors derive for the classic Kuramoto-Sivashinsky equation (e = 0) a 
second order modulation equation in the phase V' of the form ^ii — a^^x = biptxx, which is 
indeed an alternative formulation of the linearized Whitham equations; more generally, the 
Whitham equations for (1.1) consist of a system of two first-order conservation laws [NR2]. 
The signs of the coefficients a and b depend implicitly on the underlying wave train, and 
determine the spectral stability in the low frequency/low wave number regime; see again 
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[NR2]. The fact that Whitham's equations determines (a part of) the spectral stabihty 
of periodic wave trains in the low frequency regime has also been established for viscous 
conservation laws in [Se], for generalized KdV equations in [JZ4] and [BrJZ], and for the 
shallow water equations with a source term in [NRl]. 

Here, in terms of the modulation equations we have two requirements for low- frequency 
modulational stability: reality of the characteristics of the first order part (in [FST], pos- 
itivity of sound speed squared; more generally, positivity of a certain discriminant); and 
positivity of the various diffusion parameters associated with different characteristic modes 
(which are equal in the classic case [FST], but in general distinct [NRl, NR2]). 

We point out that the justification of such Whitham equations has become an important 
direction of research in the last decade, with important and fundamental results being given, 
for example, in the context of the reaction diffusion equations [DSSS] and the Shallow water 
equations [NRl] . In this paper, we carry out numerical time evolution studies that illustrate 
this correspondence with the Whitham equations; see section 2.3. In the forthcoming paper 
[JNRZl], the authors provide a rigorous justification of Whitham's equations in the context 
of viscous conservations laws, which applies in particular to equation (1.1). The approach 
proposed there extends readily to other second and higher order parabolic equations such as 
the Cahn-Hilliard equation, general fourth order thin film models, and to general systems 
of quasilinear 2r-parabolic equations such as the Swift Hohenberg equation. 

The paper is organized as follows: in the current section, we introduce the assumptions 
of our main theorem, both on the set of periodic traveling wave solutions of (1.1) and 
on the spectrum of the linearization of (1.1) about such a wave. After stating our main 
theorem we continue with section 2 where we, among other things, check numerically that 
the spectral stability assumptions of our main theorem are satisfied for some families of 
periodic traveling wave solutions of (1.1). In particular, we show that there exist bands in 
parameter space of spectrally stable periodic traveling waves and verify that these bands 
agree in particular asymptotic limits with previous numerical studies [FST, CDK, BaN]. 
Section 3 is dedicated to the proof of the main theorem of our paper: that is, that spectral 
modulational stability implies linear and nonlinear modulational stability to small localized 
perturbations. The proof is an adaptation of the analysis recently given for the St. Venant 
equations and the case of general, second order, viscous conservation laws; see [JZN, JZl]. 

For completeness, we present in Appendix A a survey of the existence theory in the peri- 
odic context in a variety of asymptotic regimes; small amplitude, near the classic Kuramoto- 
Sivashinsky equation (e— s-O, 7 = (5 = 1), and near the KdV equation (e = 1, 7 = 5 — )• 0). 
In Appendix B, we prove a nonlinear stability result, analogous to our main theorem, in 
the case of the Swift-Hohenberg equation. The nonconservative nature of this equation 
makes the nonlinear stability analysis considerably simpler than that presented in Section 
3, and the reader unfamiliar with the techniques of Section 3 may find it helpful to read the 
Swift-Hohenberg analysis as a precursor to Section 3. Also in Appendix B we conduct a 
numerical study similar to but less detailed than that given in Section 2. In particular, this 
numerical study complements and extends the analytical results of [Ml, Sc]. In Appendices 
C and D, we describe, respectively, computational statistics and the shooting algorithm 
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[BJNRZl] used to estimate the Evans function at a given frequency pair (A,^). Finally, 
in Section E, we estimate the rate of growth of perturbations of unstable waves lying near 
stability boundaries of hyperbolic vs. diffusive type. 

1.1 Equations and assumptions 

Consider the generalized Kuramoto-Sivashinsky (gKS) equation given in (1.1), written here 
in the form 

(1.2) dtu + dtu + ed^u + 6dlu + d^f{u) = 

with / G C'^{M.) and 6 > 0. The main goal of this paper is to establish that spectrally stable 
(in a sense made precise in (D1)-(D3), Section 1.1.3) periodic traveling wave solutions 
to (1.2) are nonlinearly stable with respect to small localized perturbations. As (1.2) is 
conservative, the stability analysis is an adaptation of that of [JZl, JZN] in the second- 
order viscous conservation law and hyperbolic-parabolic system cases, respectively. The 
main new aspect here is to show that the principles contained in these previously considered 
cases extend to equations with higher-order derivatives. 

It should be noted that this analysis applies also, with slight modifications, to more 
general quasilinear equations (see [JZN]), to the Cahn-Hilliard equation, and to other fourth- 
order models for thin film flow as discussed, for example, in [LP, BMS, Ho].^ Indeed, the 
argument and results extend to arbitrary 2r-order parabolic systems, so are essentially 
completely general for the diffusive case. (As already seen in [JZl], they can apply also to 
mixed-order and relaxation type systems in some cases as well.) 

We emphasize here as in the introduction above that if one considers rather the gener- 
alized Swift-Hohenberg equation 

(1.3) dtu + d^u + ed^u + 5dlu + f{u) = 0, 

considered here as a non-conservative counterpart of (1.2), the verification of nonlinear sta- 
bility of spectrally stable periodic traveling waves is considerably less complicated. Indeed, 
this is reminiscent of the distinction between general second-order parabolic conservation 
laws and systems of reaction diffusion equations; see [JZl] and [JZ2]. To aid in the un- 
derstanding of the forthcoming analysis of the conservative equation (1.2), therefore, we 
present in Appendix B the proof of the analogous nonlinear stability result for periodic 
traveling wave solutions of (1.3). Our approach there recovers and slightly extends the re- 
sults of Schneider [Sc], yielding the same heat kernel rate of decay with respect to localized 
perturbations while removing the assumption that nearby periodic waves have constant 
speed. 



^In these cases, which concerned scalar equations, periodic waves were shown to be unstable in a wide 
variety of circumstances. 
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1.1.1 Traveling- wave equation and structural assumptions 

Throughout, we consider travehng wave solutions of (1.2) or, equivalently, stationary solu- 
tions of the traveling gKS equation 

(1.4) dtu — cdxU + d^u + ed^u + Sd^u + dxf{u) = 

for some wave speed c G M. Such profiles are necessarily solutions of the traveling wave 
ODE 

(1.5) -cu' + u"" + eu'" + Su" + fiu)' = 0, 
which, integrated once, reads 

(1.6) —cu + u" + eu" + 6u' + f{u) = q, 

where g G R is a constant of integration. It follows that periodic solutions of (1.2) corre- 
spond to values {X,c,q,b) G M®, where X, c, and q denote period, speed, and constant of 
integration, and b = (6i, 62, ^3) denotes the values of (n, u' , u") at x = 0, such that the values 
of (n, u\ u") at X = X of the solution of (1.5) are equal to the initial values (61, 62) ^3)- 
Following [Se, JZl], we make the following technical assumptions: 

(HI) f e C^+^{R), K>5. 

(H2) The map H : M^ ^ R^ taking {X,b,c,q) ^ {u,u',u"){X;b,c,q) - b is fuh rank at 
{X,b,c,q), where {u,u',u"){-;b,c,q) is the solution of (1.6) such that 

{u, u , u"){0; b, c, q) = b. 

By the Implicit Function Theorem, conditions (H1)-(H2) imply that the set of periodic solu- 
tions in the vicinity of the X-periodic function u = u{-; b, c, q) forms a smooth 3-dimensional 
manifold 

(1.7) {{x,t)^U{x-a-c{^)t;p)\ {a,p)eRxU}, with ^/ C R^. 

Remark 1.1. As noted in [JZN], the transversality assumption (H2) is necessary for our 
notion of spectral stability; see (D3) in Section 1.1.3 below. Hence, there is no loss of 
generality in making this assumption. 

1.1.2 Linearized equations and the Bloch transform 

To begin our stability analysis we consider the linearization of (1.4) about a fixed periodic 
traveling wave solution u = [/(•; /3), where we assume without loss of generality X = 1. To 
this end, consider a nearby solution of (1.4) of the form 

u{x, t) = u{x) + v{x, t) 
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where v{-,t) G L'^(K) is small in some topology (to be defined precisely later), corresponding 
to a small localized perturbation of the underlying solution u. Directly substituting this 
into (1.4) and neglecting quadratic order terms in v leads us to the linearized equation 

(1.8) dtv = Lv := dx ((c — a)v) — d^v — ed^v — dd^v, a := f'{u), c := c. 

Seeking separated solutions of the form v{x,t) = e'^^v{x), v G L'^(K) and A G C yields the 
eigenvalue problem 

(1.9) Xv = Lv 

considered on the real Hilbert space L^(M). 

As the coefficients of L are 1-periodic functions of x, Floquet theory implies that the 
spectrum of the operator L, considered here as a densely defined operator on L^(]R), is 
purely continuous and, in particular, A G cr{L) if and only if the spectral problem (1.9) has 
a bounded eigenfunction of the form 

(1.10) v{x;X,0 = e'^''w{x;X,0 

for some ^ G [— vr, tt) and w G Lpgj,([0, 1]). Following [Gl, SI, S2], we find that substituting 
the ansatz (1.10) into (1.9) leads one to consider the one-parameter family of Bloch operators 

{-^«}5e[-7r,7r) acting on Ll^rii^, 1]) via 

(1.11) (L^w) (x) := e-*«^L \e'^-w{-)\ (x). 

Since the Bloch operators have compactly embedded domains Hp^j.{[0,l]) in Lp(,j.([0, 1]), 
their spectrum consists entirely of discrete eigenvalues which, furthermore, depend contin- 
uously on the Bloch parameter ^. It follows then by standard considerations that 

0-L2(M)(-^) = U CrL2„([0,l]) (^«) ; 

Ce[-7r,7r) 

see [Gl] for details. More precisely, the spectrum of L can be characterized by the union of 
countably many surfaces A(^) such that A(^) G 0"(L^) for each ^ G [— vr, tt). 
Given a function u G L^(M) we now recall its inverse Bloch-representation 

(1.12) u{x)= r e'^-''u{i,x)di, 



where n(^, x) := ^^ e^'^*'^^u(^ + 2TTk) are 1-periodic functions and u{-) denotes the Fourier 
transform of u, defined here as 

2tT ./in. 
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Denote hy B : u >-^ ii the Bloch transform operator taking u G L'^{]^) to its Bloch transform 
u £ L^([— vr, vr); Lpgj.([0, 1])).^ For a given linear operator L with 1-periodic coefficients, 
it can readily be verified that B (Lu) {^,x) = L^[u{^,-)]{x) hence the associated Bloch 
operators L^ may be viewed as operator-valued symbols under B acting on Lpgj.([0, 1]). 
Furthermore, the identity B (e^^u) (^, x) = e^^^ [^(C) ")] {^) naturally yields the inverse rep- 
resentation formula 

(1.13) e^*no= re^«-^e^«M^o(e,-)](x)dC 

J — IT 

relating behavior of the linearized equation (1.9) to that of the diagonal Bloch operators 

The representation formula (1.13) suggests that bounds on the linear solution operator 
e acting on some function in L^(]R) may be obtained from bounds on the associated 
operators e^^* acting on Lpgj.([0, 1]). To facilitate these estimates, we notice by Parseval's 
identity that the Bloch transform u i— )■ B{u) = u is an isometry in L?, i.e. 

/TT p\ 
/ \u{i, x)\^ dx dC = Mu\\l2^[_^^^),L2([o^i])), 
-TV Jo 

More generally, we have for any one-parameter family of 1-periodic functions /(^, •), ^ G 
[— TT, tt), the generalized Hausdorff- Young inequality 

^ ^^■^^^•^^^",,,^ <(2-)^II/IIm([-.,.),l.([o,i])) 



(1.15) / e*^7(e,-)rfe 



— 7r 



valid for any q < 2 < p and ^ + ^ = 1; this is readily obtained interpolating Parseval's 
identity with Triangle Inequality (case q = 1, p = oo), see [JZl]. It is in the context of the 
above functional framework that we will obtain our crucial linearized estimates; see Section 
3.2 below. 

1.1.3 Spectral stability assumptions 

Taking variations in a neighborhood of u along the 3-dimensional periodic solution manifold 
(1.7) in directions for which the period does not change, we find that the generalized kernel 
of the Bloch operator Lq is at least two-dimensional by assumption (112); see [NR2] for 
more details. Following [JZl], then, we assume along with (II1)-(II2) the following spectral 
diffusive stability conditions: 

(Dl) a{L) C { A G C I Re(A) < } U {0}. 

(D2) a{L^) C { A G C I Re(A) < -6*1^12 }, for some 6^ > and any ^ G [-7r,7r). 

(D3) A = is an eigenvalue of Lq of algebraic multiplicity two. 



^Here, and elsewhere, we are adopting the notation ||/||L5([-7r,7r),i,p{[(),i])) := (fl^ Wfi^^ ■)\\lp ao.i])'^y 



i/g 



1 INTRODUCTION 11 

Assumptions (H1)-(H2) and (D3) imply that there exist two eigenvalues of the form 

(1.16) A,-(0 = -m,-e + o(|e|) 

of L^ bifurcating from A = at ^ = 0; see Lemma 3.1 below. Following [JZl], we make the 
following non-degeneracy hypothesis: 

(H3) The coefficients Uj in (1.16) are distinct. 

Hypothesis (H3) ensures the analyticity of the functions Aj(-); again, see Lemma 3.1 below 
for details. 

2 

Continuing, we notice that when f{u) = ^ the Galilean invariance of (1.1), along with 
assumptions (H1)-(H2) and (D3), implies that the zero-eigenspace of Lq admits a non-trivial 
Jordan chain of height two [NR2]. Indeed, in this case it is readily checked that the map 

Qcu{x, t) = u{x — ct, t) + c 

sends solutions of (1.4) to solutions of (1.1) so that, in particular, variations in wave speed 
are period-preserving. It follows that the periodic solution manifold (1.7) can be parame- 
terized as 

(1.17) { (x,i) h^C/(x-a-ct;c,6) I (a,c,/3) gR2 xX }, with X C M. 

that is, wave speed serves as a non-degenerate local coordinate on (1.7). Noting that 
{dt — {L — cdx))Qc = Qc{dt ~ L), we remark that variations along the periodic manifold 

(1.17) correspond to solutions of the (formal) linearized equation (dt — {L — cdx))v = and 
find, setting u = [/(•; c, 6), that 

{dt -{L- cdx)) {dclGcUi-, c, 6)]|,=5) = 0. 

Since 5c[^cii(-; c,6)]|c=c = Qc{-tUx) + ^3cC/|(c,6)=(c,fe) and 9ct/|(cf,)^(g5) is periodic of period 
1, this verifies the existence of a Jordan chain of height two (recall, the height is at most 
by (D3)) ascending above the translation mode Ux- 

For general nonlinearities /, however, the existence of such a Jordan block does not 
immediately follow from assumptions (II1)-(H2) and (D3). Indeed, setting /3 = {fii, P2) £ '^ 
in (1.7), we find, taking variations along (1.7) in I3j near u = U{-](5), that 

(1.18) LdpV\p^-p = -dp.c\p^-pUx. 

As the function df^j^X{^)df^2^{-; (3) — dii^X{J3)dii^U{-]J3) is periodic of period X{J3) = 1^, we 
see that the existence of a non-trivial Jordan block requires the additional non-degeneracy 
hypothesis^ 

^Here, we notice by (D3) that P can not be a critical point of the period X. Indeed, if dX(j3) = then 
dpiU{-;P) and dp-^U{-;l5) would both be 1-periodic, hence A = would be an eigenvalue of Lq of algebraic 
multiplicity three by (1.18). 

"^For more general systems, this requirement would be replaced with the condition that the gradients of 
the wave speed and period be linearly independent at (3. 
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For a much more thorough and precise discussion on this topic in the context of general 
systems of conservation laws, please see the forthcoming work [JNRZl]; see also Remark 1.3 
and Lemma 3.1 below. Throughout this work, we assume that such a Jordan block exists 
by enforcing our final technical assumption: 



(H4) The eigenvalue of Lq is nonsemisimple, i.e., dimkerLg = 1 



5 



Remark 1.2. The coefficients Oj are the characteristics of an associated Whitham averaged 
system formally governing slowly modulated solutions; see [NR2]. Thus, as discussed in 
[0Z4, JZl, JZN], (Dl) implies weak hyperbolicity of the Whitham averaged system (reality 
of Oj), sometimes used as a formal definition of modulational stability; condition (H4) holds 
generically, and corresponds to the assumption that speed c is non-stationary along the 
manifold of nearby stationary solutions, see Lemma 3.1; and condition (D2) corresponds to 
"diffusivity" of the large-time (~ small frequency) behavior of the linearized system, and 
holds generically given (H1)~(H4), (Dl), and (D3), see also [SI, S2]. 

Remark 1.3. As noted in [NR2, JNRZl], the conservative structure of (1.1) implies that, 
assuming (HI) and (D3), one can take /3 = {X,M) in (1.7), where M = X~^ L U{s;f3)ds 
denotes the mean mass of the associated wave. Indeed, following the computations of [0Z3] 
we find under these assumptions that for |A| ^ 1 

where F 7^ is a constant, k = X^^ denotes wave number and Dq{-) represents the Evans 
function^ for the Bloch operator Lq, from which our claim follows by assumption (D3). 
Using this parameterization of the manifold (1.7), condition (1.19) reduces to 

corresponding to linear phase coupling in the language of [JNRZl]. In particular, notice 

2 

that in the case f{u) = \ o'^^ can use Galilean invariance to show 

yielding another verification of (H4), assuming (H1)-(H2) and (D1)-(D3), in this case. 

1.2 Main result 

With the spectral assumptions discussed in the previous section, we are now prepared to 
state the main theorem of this paper. 



"'The degenerate case where (1.19), hence (H4), fails can be treated as in [JZ3]. 

^Recall that the algebraic multiplicity of the root of the Evans function at A = agrees with the algebraic 
multiplicity of the eigenvalue A = of the associated linear operator; see [Gl]. 
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Theorem 1.1. Let u be any steady 1-periodic solution of (1.4) such that (H1)-(H4) and 
(D1)-(D3) hold. Then there exist constants Eq > and C > such that for any uq with 
11^0 — ^llLi(M)nH^(R) ^ ^0) where K is as in assumption (HI), there exist u a solution of 
(1.4) satisfying u(-,0) = uq and a function tp{-,t) G W^'°°(R) such that for all t > and 
2 < p < oo we have the estimates 



--(1--) 

\\u{-+^i-,t),t)-u\\LP(R)<C{l + t) H ^^ ||M(-,0)-u(-)||Ll(IR;)nH^(R) 

(1_20) ||(V't,V'x)(-,i)llLP(R)<C(l + t) H -n\u{;0)-u{-)\\L^R)nH^iR) 

M- + ^(•, t),t) - u\\hk^^) < C (1 + t)-3 ||n(., 0) - ^(OILi 
'^t, i^.){; i)llH^(R) < C (1 + t)-3 ||n(., 0) - u{.)\\^. 



Moreover, we have the L-'^(M) n H^(K) — )■ L°^(R) nonlinear stability estimate 

(1-21) ll^(')*) - ^IIl°°{R) ' II^(">OIIl°o(R) ^ C" ||'u(-,t) - u\\L^(R)nH^{M.) 

valid for all t > 0. 

Theorem 1.1 asserts asymptotic L^ n H — )■ L^ convergence of the modulated solution 
u{- + tlj{-,t),t) toward u, but only bounded L^ n H^ — t- L°° nonlinear stability of the 
unmodulated solution u(-) about u. 

Remark 1.4. It may seem more natural to introduce ip via v{x, t) = u{x, t) — u{x + ^l^{x, t)) 
However, in doing so one introduces in the equation for v terms involving only tp and thus 
not decaying in time. For this reason we work instead with v{x, t) = u{x + ijj{x, t),t) — u{x), 
that is, u{x, t) = u{Y{x, t)) + v{Y{x, t),t), where Y is such that Y{x, t) + tp{Y{x, t), t) = x, 
Y{y + ilj{y,t),t) = y. Notice that we insure the existence of such a map Y by keeping, for 
any t, \\ip{-,t)\\icxi(^-^ bounded and ||V'x(-,i)||L°°(R)- It should be stressed, however, that 

y(x,t)=X-V(x,t) + 0(||^(-,t)||L-(M)||V':.(-,t)|U-(M)) 

SO that we are not so far from the natural (but inappropriate) approach; see [JNRZl] for a 
detailed discussion. Notice, moreover, that introducing the map Y above enables one to go 
back to the original unknown u{x, t). 

The proof of Theorem 1.1 is presented in Section 3. As noted earlier, the proof of 
the corresponding theorem for periodic traveling wave solutions of the (non-conservative) 
Swift-Hohenberg equation (1.3) is considerably less complicated than that of Theorem 1.1; 
the reader unfamiliar with the methods of this paper may thus wish to consult Appendix 
B before proceeding to the proof of Theorem 1.1. 

While it is generally accepted that the structural hypotheses (H3)-(H4) and (D3) should 
generically hold, it is unclear a priori whether there exist periodic traveling wave solutions of 
(1.2) that satisfy the spectral diffusive stability assumptions (D1)-(D2). While an analytical 
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verification of these conditions seems beyond the scope of our current machinery, due to the 
complexity of the governing equations, we provide a numerical study of this topic in Section 
2 below giving convincing evidence of existence of stable bands of periodic waves. We view 
the rigorous numerical or analytical proof of spectral stability in some interesting regime as 
the primary remaining open problem in this area of investigation. See in particular [BaN] 
and [JNRZ4] for progress toward a numerical proof in the KdV limit |(7, (5)| — s- 0. 

2 A numerical study: spectral stability and time evolution 

In this section, we address numerically the issue of spectral stability of periodic traveling 
wave solutions of (1.1). Our numerical study is based on complementary tools, namely. 
Hill's method and Evans function computations; see [BZ, BJNRZl]. On the one hand, 
we use SpectrUW numerical software [CDKK], based on Hill's method, which is a Galerkin 
method of approximation. More precisely, the coefficients of the linear operators L^ (defined 
by (1.11)) together with periodic eigenvectors are expanded by using Fourier series and then 
a frequency cutoff is used to reduce the problem to finding eigenvalues of a finite matrix. It 
is known that Hill's method converges "globally" very fast, but with rates not particularly 
guaranteed in fixed areas; this method is used here to obtain a good localization of the 
spectrum. In order to test rigorously the spectral stability, we use an Evans function 
approach: after bounding the region where unstable spectra might exist, we compute a 
winding number along various contours to localize the spectrum. In some sense, Hill's 
method helps to find suitable contours to test the spectral stability. To analyze stability 
near the origin, we also compute a Taylor expansion of the critical eigenvalues bifurcating 
from the origin at the ^ = state. The method of moments (to be described in section 
2.1.5) then guarantees the stability of these critical eigenvalues in the neighborhood of 
(A,.^) = (0,0). While most of the results of this section are numerical, we point out that 
the numerics are well-conditioned with rigorous error estimates; see [BJZ, Zl, BHZ2] for 
more details. 

Throughout this section, both for definiteness and in order to compare with some pre- 
vious numerical results (see e.g. [BaN, EST, CDK]), we consider (1.1) unless otherwise 

2 

specified with the specific nonlinearity f{u) = A^, A > 0, i.e. the equation 

(2.1) dtu + 6d'^u + ed^u + ^d^u + AudxU = 

with various values of e, 6, 7, A. If one introduces the characteristic amplitude U, length 
scale L and time scale T = L/AU, equation (2.1) is alternatively written as 

(2.2) diu + ud,u + j^dlu + j^dlu + j^d^u = 0. 

As a result, one can always assume A = 1 in (2.1). Depending on applications and numerical 
purposes, one can remove two of the three other parameters in the study of (2.1). We have 
chosen here two reduced parameterizations of the problem. On the one hand, we first choose 
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some e then, setting e = e/AUL'^ and 7 = KUL^ into (2.2), one finds that it is equivalent to 
choose 7 = 1 and e is a fixed parameter and let 5 vary as a free parameter in (2.1). On the 
other hand, letting L^ = ^/5 and U"^ = (e/AL^)^ + ((5/AL)^, one finds that it is equivalent 
to choose^ 7 = (5 and 5^ + e^ = 1 in (2.1). This latter case is particularly of interest since it 
is found in thin film theory^ and since all equations (2.1) can be reduced to this particular 
form through rescaling; that is we cover all the values of parameters e, (5, 7. However this 
form of the equation may lead to a singularly perturbed problem as (5 — )■ 0; that is why we 
have also focused on the first form of (2.1) with 7 = 1 even if, to cover all cases, we need 
then to introduce two families, one with e fixed to and the other one with e fixed to some 
arbitrary positive value. 

The plan of the study is as follows. As a first (testing) step, we consider the spectral 
stability case 7 = 1 and e fixed. In this test case, we give all the details of our numerical 
approach: we first compute a priori bounds on possible unstable eigenvalues and then 
describe our numerical Galerkin/Evans functions based approach. In order to compare our 
results with the existing literature, specifically the results in [FST], we first carry out a 
numerical study in the case of the classical Kuramoto-Sivashinsky equation (e = 0, 7 = 1). 
We follow this with a numerical study for the generalized KS equation when e = 0.2 and 
7 = 1. Next, we switch to the "thin film scaling" 7 = (^, 5^ + e^ = 1 and show numerically 
that there are (sometimes several) bands of stable periodic traveling waves for all choices 
of the model parameters, comparing to previous numerical studies of [BaN] and [CDK] 

We end this section with a discussion of time-evolution studies nearby various periodic 
traveling wave profiles numerically computed and compare with the dynamics predicted by 
the associated Whitham averaged system. 

2.1 Spectral stability analysis: 7 = 1, e fixed 

2.1.1 Continuation of periodic traveling waves 

Let us consider periodic traveling wave solutions of (2.1) in the simple case 7 = 1 and e 
is fixed. After integrating once, traveling waves of the form u{x^ t) = u{x — ct) are seen to 
satisfy the first order nonlinear system of ODE's 

u \' ( u' 

(2.3) I n' = u" 

u" I \ cu- eu" -5u' -y^ + i 

for some constant of integration (7 G M. Noting that we can always take c = by Galilean 
invariance, we find two fixed points in the associated three-dimensional phase space given 



'^Notice that, applying the transformation {x, u) 1— > (— x, —u) if necessary, we may always take e > in 
(2.1). Thus, here we set e = ^/l-S"^ or <5 = Vl - £^- 

®Where the equation is designed to investigate formation of patterns of size 1 at a threshold of instability 
and therefore is scaled as 7 = 5 so that, for constant states, the most linearly- unstable perturbation has 
frequency 1. 
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by 

U-iq)= \ I and U+{q) 





Following the computations of Section A.l, we find that the fixed point Uj^{q) undergoes 
a Hopf bifurcation when ei^HopfCfz) = \/2^j see Figure 1(a) when e = 0.2. Moreover, near 
this bifurcation point we find a two dimensional manifold of small amplitude periodic orbits 
of (2.3) parameterized by period and the integration constant q. Observe that due to the 
presence of the destabilizing second order and stabilizing fourth order terms in (1.2), an 
easy calculation shows that the equilibrium solution U+{q) = \/2q^ and in fact all constant 
solutions, are linearly unstable to perturbations of the form v{x,t) = e''**+*^^, A G C, A; G M, 
with \k\ <C 1, while linearly stable to such perturbations when |A;| ^ 1. As a result we 
find that all small amplitude periodic orbits of (2.3) emerging from the Hopf bifurcation 
correspond to spectrally unstable periodic traveling waves of (2.1). 

Not surprisingly, we are able to numerically continue this two parameter family of small 
amplitude periodic orbits to obtain periodic orbits of (2.1) with large amplitude. Further- 
more, in the natural three-dimensional phase space of (2.3) and for e fixed, we are able to 
find a one-parameter family of periodic orbits of (2.3) with fixed period X by continuing 
the periodic orbits emerging from the Hopf point Uj^{q) with respect to the integration 
constant g; see Figure 1(b) for e = 0.2 and X = 6.3. Notice that in order to keep the period 
X fixed, one naturally has to vary the modeling parameter 5 with respect to q. 

What is perhaps surprising, however, even though by now well known [FST, CDK], is 
that there appears numerically to be a small band of spectrally stable periodic traveling 
waves within this one-parameter family, outside of which all waves are spectrally unstable; 
see Figure 2 and Figure 3. Our immediate goal is to substantiate this claim with careful 
numerical procedures completed with rigorous error bounds. Although, due to machine 
error, the justifications we discuss do not constitute a proof, given the rigorous error bounds 
involved, our claims could in principle be translated into a numerical proof with the help 
of interval arithmetic. This would be an interesting direction for further investigation. 

Remark 2.1. As well explored in the literature (see for example [CDK]), there occur a 
number of period doubling bifurcations as parameters are varied, which might in principle 
be difficult to follow by continuation. However, this difficulty is easily avoided by rescaling 
so as to fix the period, while letting other parameters vary instead. This is our reason for 
computing with fixed period as we do. 



2.1.2 High-frequency analysis 

We begin our study of the spectrum of the linear operator L, obtained by linearizing (2.1) 
about a fixed periodic traveling wave solution, by eliminating the possibility that the as- 
sociated Bloch operators L^ admit arbitrarily large unstable eigenvalues. As the resulting 
estimates are independent of the Bloch frequency ^, this allows us to reduce our search for 
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(a) 




(b) 




Figure 1: (a) Here, we plot a few of the trajectories in the phase portrait of the (2.3) for 
£ = 0.2, q = 0.04, and (5 ~ 1. The periodic trajectory bifurcating from the nearby unstable 
equilibrium solution U+{q) is plotted in bold, and has period X = 6.3. (b) This figure 
depicts the one-parameter family (up to translations) of periodic orbits of period X = 6.3 
and wave speed c = emerging from the unstable equilibrium solution U^{q). Here, e = 0.2 
and q (and accordingly 5 = 6{q)) varies from 1 (corresponding to the orbits near the Hopf 
bifurcation) to 30 by unit steps. The bold periodic orbits near the center correspond to the 
lower and upper stability boundaries: in particular, all periodic traveling wave solutions 
of (2.1) corresponding to periodic orbits of (2.3) between these bold orbits are spectrally 
stable with respect to small localized perturbations. 
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Figure 2: A numerical sampling of the continuous spectrum of the linearization of (2.1) 
when e = 0.2, plotted here as Re(A) vs. Im(A), about the periodic traveling wave profiles 
on the periodic manifold depicted in Figure 1(b) for (a) q = 1, (h) q = 4, (c) q = 5, (d) 
q = 6, (e) q = 7, (f) q = 12, (g) q = 13.8, and (h) q = 30. Moreover in (i) we zoom in 
near the origin in the spectral plot given in (c). All pictures here were generated with the 
Galerkin based SpectrUW package developed at the University of Washington [CDKK]. 
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Figure 3: Here, we zoom in near the origin in the A-plane near, but just below, the lower 
stability boundary when considering 6.3-periodic stationary solutions of (2.1) in the case 
e = 0.2. The three figures correspond to (a) q = 4.4, (b) q = 4.5, and (c) q = 4.6. We 
find that prior to stability, a real anti-periodic eigenvalue (corresponding to ^ = — vr, and 
marked with a large star) crosses through the origin and then the rest of the unstable 
spectrum crosses the imaginary axis away from the origin. This crossing of a real anti- 
periodic eigenvalue corresponds to a period doubling bifurcation; see [KE] for more details 
on period- multiplying bifurcations for KS. The spectral curves were generated by computing 
the roots of the Evans function using the method of moments described in Section 2.1.5. In 
particular, the transition to stability is not signaled by an eigenvalue crossing through the 
origin. 
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unstable spectrum of L into a fixed compact domain in A. For generality, here we consider 
(1.1) with a general (smooth, say) nonlinearity / rather than the quadratic nonlinearity 
leading to the formulation in (2.1) and 7 = 1. 

Letting -u be a fixed X-periodic solution of (1.2), for each ^ € [—tt/XjIt/X) we must 
study the eigenvalue problem 

An + {d^ + iO {f'{u)u) +eid^ + i^f n + <5 (9^ + i^f u + {d^ + i^^ u = 

considered on L'^^^{[0,X]) We restrict ourselves to rational Bloch frequencies, that is to 
cases where S,X /{2tt) £ Q, (recall that the spectra of L^ is continuous in ^) and thus this is 
equivalent with studying the family of eigenvalue problems 

(2.4) Xu + {f'{u)uy + eu'" + 6u" + u"" = 

considered on Lpgj,([0, nX]) for each n G N. In this latter formulation, we have the following 
estimates on the modulus of possible unstable eigenvalues. 

Lemma 2.2. Let n G N 6e arbitrary, and suppose that there exists a function u G L?^^^{[0, nX]) 
which satisfies (2.4) for some A G C. Then A must satisfy the estimates 

(2.5) Re(A)<^||r(n)n'|Lo.([o,x]) + J, 
and 



(2.6) ReA + |ImA| < ^ (^\\f"iu)u'ho.^[oM + \\f'imU(^io,x]) + ^ 



2+(l + 2£T 



2 

In particular, the above estimates are independent of n. 

Proof. By taking the real part of the complex scalar product of (2.4) with u, one obtains 

1 /"nX 

(2.7) M>^)Ml2^^^[o,nX]) = "2 J^ (/'(^))>l'(2^)^2; + 5||n'||22^^(p^„^]) - lk"|li2^^([o,„x])- 
Then, use the Sobolev interpolation bound 

with (7 = 2 to deduce (2.5). Next, by taking the imaginary part of the complex scalar 
product of (2.4) with n, one obtains 

then 

\^'^W\\Ml2^^(^lO,nX]) ^ ll/'(^)llL-([0,X])hllL2,,([0,nX])lk'llL2,,{[0,nX]) 

+ 4u\\Lj,,{[0,nX])h"\\Ll,A[0,nX])- 
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By using Young's inequality, one finds 

valid for any D > 0. By using again the above Sobolev interpolation bound, one finds 



(2.8) 



1 + eD e \ „ 2 



where now C,D > are arbitrary constants. From (2.7), one also deduces that 

(2.9) M^)\Ml2^^^[o,nX]) ^ 2^\\f"i^)^'\\L^{[0,X]) + ^^)h\\LH[0,nX]) " 2^\^"\\lmo,nX])- 



Choosing now D = 2e, C = I + De and adding (2.8) and (2.9) yields 



Re(A) + |Im(A)| < ^(||/"(t2)n'||i^([o,x]) + ll/'(^) 



,2 ,2 (l + 2e2)2 



lL°°([0,X]) "^'^ "T 2 

as claimed. D 

Remark 2.3. Notice that by choosing different constants C and D above, along with a 
different choice of the Young's scaling in (2.9), various other bounds of the form (2.6) can 
be obtained which could be beneficial in certain parameter regimes. Specifically, given any 
constants Ci, C2, C3, a, /3,A > such that 

^ 11/1 /^^2|| W/-NII , e/3C'3 e 

2cr-l + ^(4a"^^")"^-M) + ^ + 2^ 

we obtain the inequality 

Re(A)+^|Im(A)| < l||/"(n)^'||^.([o,^]) + ^ + | {\\f'{mL^([o,x]) (" + ^) + ^) • 
The inequality (2.6) corresponds to the choices 

C, = 6, C2 = C3 = ^^, a=\\f'{u)\\L^iio,x])' '^ = 2e' ^ = 1' 

while an example of an alternative bound, which may be useful for large e > 0, is given by 
(2.10) 



Re(A) + ^|Im(A)| < ^ (|ir(n)n'|L.([o,x]) + <^' + ^ (^ll/'(^)llio.([o,x]) 
corresponding to the choices 



+ 3? 



3^2 -^ 2 1 

Ci = 6, C2 = C3 = —, a= -||/'(^)||i.c([o,x])' ^=£' ^"2e2- 
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By Lemma 2.2, for fixed modeling parameters e and 6 and a periodic profile u of (1.2), 
there exists a real number Rq = Ro{e,6,u) such that the unstable (and marginally stable) 
spectrum of the associated linearized operator L satisfies 

(2.11) cr (L) n { A G C I Re(A) > } C B{0, Rq). 

This reduces to a compact set the region in the spectral plane where we need to search for 
unstable spectra of L. In order to verify the spectral stability hypotheses (D1)-(D3) and 
(H3) we use the periodic Evans function and basic analytic function theory, the necessary 
details of which we review in the next section. 

2.1.3 Low- and mid-frequency analysis: Evans function and numerical methods 

We now begin our search for unstable spectra of L in the ball B{0, Rq), where Rq is given by 
the high-frequency bounds in the previous section. To this end, we use a complex-analytic 
function, known as the periodic Evans function [Gl], that is well-suited to the task at hand. 
First, let n be a fixed X-periodic traveling wave solution of (2.1), X arbitrary, and notice 
that the associated spectral problem (1.9) can be rewritten as a first order system of the 
form 

(2.12) Y'(x;A) =M(x,A)Y(x;A), Y G C^ 

and that A G o"(L) if and only if (2.12) admits a non-trivial solution satisfying 

Y{x + X;X) = e'^^Y{x;X), Vx G M 

for some ^ G [—tt/XjTt/X).^ Letting ^(x,A) be a matrix solution of (2.12) with initial 
condition ^(0, A) = I for all A G C, we follow Gardner [Gl] and define the periodic Evans 
function for our problem to be 

(2.13) D{X,0-=det('^{X,X)-e'^^l), {X,O^C^, 



where for later convenience we allow for possibly complex Bloch frequencies ^ in the above 
definition. By construction, then, A G o"(L) if and only if there exists a ^ G [—ir/X, tt/X) C 
M such that D{X,^) = 0. More precisely, for a fixed ^ G [—tt/X,tt/X) the roots of the 
function D(-,^) agree in location and algebraic multiplicity with the eigenvalues of the 
Bloch operator L^ [Gl]. 

Since ]HI(-,A) clearly depends analytically on A, it follows that the Evans function is a 
complex-analytic function of both A and S^, hence our search for unstable spectra of L may 



®Here, we have allowed the domain of ^ to depend on the period X. This choice is to ensure that, when 
linearizing about a constant state, the dispersion relation D{X, £^) — agrees with that obtained by the 
standard Fourier transform, where D denotes the periodic Evans function defined below in (2.13). This 
differs from the convention of Gardner [Gl] and some others in which the Evans function is parametrized 
instead by the Floquet number i^ := e*^^ G S*^, a quantity with domain independent of X. 
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by reformulated as a problem in analytic function theory. Indeed, for a fixed ^ the number 
of eigenvalues of L^ within the bounded component of a closed contour F, along which 
D{-,C)\r is non- vanishing, can be calculated via the winding number 

^^' ' 2m Tt 0(A,{) 

Given a fixed contour T then, notice as n(-; T) is a continuous function with integer values, it 
is locally constant and hence can only change values when a root of D{-,^) crosses through 
r. With these preparations in mind, we now outline the general scheme of our numerical 
procedure to verify the spectral stability hypotheses (D1)-(D3) and (H3). As noted below, 
for numerical efficiency, the following steps are sometime carried out in a slightly different 
order. Of course then, there are several consistency checks that follow to justify various 
assumptions made. For similar reasons, it is at times advantageous to execute a particular 
step in a slightly different way than suggested here; more will be said about this below. 
Nevertheless, the general idea of the numerical procedure is contained in the following 
steps, to be completed for fixed. 

Step 0: First, numerically determine set V := {q\ (?|Z)(0, 0) = 0} by looking for sign changes^^ 
of 9|-D(A,0) for A > 0, where here -D(A,0) denotes the Evans function for .^ = 
obtained by linearizing about the periodic profile u with parameter q. 

Then away from this degenerate parameter set T>, we may complete the following steps. 

Step 1: First, determine using the bounds of Lemma 2.2 an i?o > satisfying (2.11). 

Step 2 : Let < tq < Rq, ^q := B{0, Rq) \ S(0, ro), and set 

ro:=a(Oon{AGC|Re(A)>0}), 
^' ' Fi := 9(5(0, i?o)n {AG C|Re(A) >0}); 

see Figure 4. Now, verify that for some < /cq ^ 1 the following two conditions hold: 

(a) n(^,Fo) = for aU ( G [-tt/X,7:/X). 

(b) n(e,Fi) = Oforfco<|e|<7r/X; 

see Figure 5 below. 
Step 3: Verify that for some tq < ri < Rq, n{^, dB{0, ri)) = 2 for |^| < /cq. 
Step 4: Finally, Taylor expand the two critical modes for |^| <^ 1 as 

(2.15) A,-(0 = «ie + /3,^' + O(|^|'), i = l,2 

and numerically verify that aj G Mi, the aj are distinct, and Re(/3j) < 0. Notice that 
such an expansion exists by Lemma 3.1. 



^Since the map q n> dxD{0, 0) is real valued, this is sufficient to find the zeros. 
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If the criteria of each of the above steps are satisfied, then we conclude (numerically) 
that the wave is stable. We emphasize here the simplicity and generality of the above 
numerical protocol, hence being directly applicable to many model problems. In particular, 
the above numerical method does not require an often lengthy and problem specific spectral 
perturbation expansion near the neutral stability mode (A, ^) = (0,0). 



(a) 




(b) 
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Figure 4: Illustrations of the curves Fq (a) and Fi (b). These pictures are not drawn to 
reflect the actual dimensions used in the study, but to be a visual aid. 
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Figure 5: Sample Evans function output, demonstrating the calculation of n(^, Fq) in Step 
2 of the above numerical protocol. We view in (a) a domain contour of form Fq defined in 
Step 2 above on which we evaluate the Evans function, and in (b) the associated range for 
^ = 0. Here, e = 0, g = 5, X = 6.3, and b = 1.4627. In particular, we see that n(0,Fo) = 
in this case. 



Remark 2.4. The set T> introduced in Step above is precisely the parameter set in which 
hypothesis (H2) is valid; see [NR2, JNRZl] for details. In particular. Step is a useful test 
for co-periodic instability corresponding to bifurcation of periodic solution to a new branch. 
Prom a numerical standpoint, it should be noted that Step 3 of the above procedure will 
fail, producing an infinite loop in the algorithm, on the set D. In general, given a numerical 
tolerance < 7 <C 1, the algorithm is only guaranteed to be well-conditioned off of a 
7-neighborhood of T) since this guarantees we are away from a degenerate boundary case 



2 A NUMERICAL STUDY: SPECTRAL STABILITY AND TIME EVOLUTION 24 



where d^D{0, 0) is too small. This will produce sharp stability boundaries in the region 
where (9^D(0,0) is small. Such degenerate boundaries, while seemingly not occurring in 
the present case, do occur in the non-conservative Swift-Hohenberg equation studied in 
Appendix B. 

Remark 2.5. It is clear that some care must be taken in choosing the low-frequency cutoff 
ko in the above procedure, as it must be chosen small enough to ensure that the conditions 
on ttj and f3j in Step 4 are sufficient to guarantee that Re(Aj(^)) < for |^| < ko; that is, 
the cubic order remainder term must be small enough for |.^| < /cq as to not dominate the 
O(I^P) terms determined by Taylor expansions (2.15). To analyze how small /cq must be, 
assume |^| < ko and note that Taylor's theorem implies 



A,(e) = a,^ + f3,e + y / (1 - sfX'fisOds. 



Letting K ^ ko and setting^^ Mj = max|^|=;^ l^j(C)l) where ko is given in Step 3 above, it 
follows by basic interior estimates that for all s G [0, 1] 



I a;" (.01 



2vrJ JdB{o,K) (C - sO" 



< 



6MjK 
[K^ko] 



4" 



In particular, upon verifying that Re(aj) = and Re(/3j) < for j = 1, 2, we can guarantee 
Re(Aj(0) < for |^| < ko provided that ko satisfies the inequality 

(2.16) ko < ^ • min ' ^^^^' 



K j=i,2 Mj 

This serves as a consistency check in the above numerical procedure: if the chosen ko does 
not satisfy the stated bound then one must choose a new ko which does satisfy the bound 
and repeat Steps 1-4 with this new ko- 

2.1.4 Numerical determination of Taylor coefficients 

Concerning Step 4 in the above procedure, we now explain how the coefficients aj and f3j in 
(2.15) can be obtained from the Evans function. First, assuming that Steps 1-3 have been 
verified, the analyticity of the Evans function implies that D{-, ■) can be Taylor expanded 
to third order as 

D{\, = aoA^ + aiA^ + aaC^ + ogA^ + a^X^^ + a^Xf + aef + 0{\X\^ + \^f) 

where the coefficients oq , . . . , oe G C can be calculated from the Evans function via the 
formula 

(2.17) dldW(0,0) = ^-^ (( i D(fi,k)n''-^k-'-^dtidk 

'^ 47r^r! JdB(o,h) JdB(o,h) 



11 
below for details 



In our numerics, we utilized the method of moments to determine the number Mj; see Section 2.1.5 
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where h > is sufficiently small. Substituting the expansions in (2.15) into the equation 
D{Xj{^),^) = and equating second-order terms, we find that the aj are roots of the 
polynomial agz'^ + aiz + 02 = 0. In particular, we can take 

(2.18) a, = -«i + (-lP^V«^- 40002 

noting that (strong) hyperbolicity of the associated Whitham averaged system is equivalent 
to the discriminant condition 

(2.19) a?-4aoa2<0, 

which is easily checked. Furthermore, equating third-order terms in D(Aj(^),.^) = implies 
that 



(2.20) (3j 



2aQaj + oi 



where we note that the denominator is non-zero by the definition of the aj and the discrim- 
inant condition (2.19). 

2.1.5 Method of moments 

Finally, we describe how to numerically obtain bounds on the numbers Mj := maxi^i^^ l'^i('^)l- 
From the above discussion, the computation of these quantities is important to determine 
the sign of eigenvalues in the neighborhood of the origin. Here, the idea is to use the Evans 
function to track the spectral curves Xj{k) as functions of k by using what we refer to as 
the method of moments. We need to locate two eigenvalues: to this end, we define the first 
and second moments of the spectral curves Xj{^) as 

^liO ■■= E HO m2{0 := E ^^•(^)'- 

j=l,2 j=l,2 

Note that these moments can be directly computed through the Evans function via the 

formula 

1 / \'dxD{\,i) 



for i? > sufficiently large so as Aj(^) G B{0,R) for all |^| = K. Once the moments mi{S,) 
are computed for a fixed ^, by using a simple Stieltjes trapezoid method, for example, the 
corresponding numbers Xj{^) can be easily recovered via the formula 



A (n = "^^^^) + (-l)V2m2(g)2 - mi (02 
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from which the numbers Mj can easily be determined. This method of moments readily 
extends to an arbitrary number of eigenvalues; see [BZ, BJNRZl]). In order to avoid the 
ambiguity for the choice of Ai and A2 above, we substitute (2.16) by 



01 

KM j=i,2 ' 



^0 < ^^77 — ■ ™ii\ \^^{Pj] 



where M = maxj=i^2 max|^|=^ l^^ilOI- 

2.1.6 Numerical results 

The above numerical procedures provide a direct and well-conditioned way to numerically 
verify the spectral stability hypotheses (H3) and (D1)-(D3). We now apply this procedure 
to two cases. 

First case: 7 = 1,£ = 

The case e = corresponds to the "classic" Kuramoto-Sivashinsky equation, which has 
received a great amount of numerical and analytical attention over the years. In particular, 
stability boundaries are well known for this equation. We state our results first in this 
special case for verification of our method against these previously known results. Let us 
now carefully describe the implementation of the numerical procedure for the special case of 
(2.1) for the periodic traveling waves parameterized by g = 5, and X = 6.3. For numerical 
efficiency, we proceed in our explicit example by slightly modifying both the content and 
order of the steps of the above numerical procedure. As a first step, we notice that when 
e = and q = 5, a periodic orbit of (2.3) with period X = 6.3 exists when 6 = 1.4627, and 
so we fix 6 at this value throughout this discussion. The corresponding solution u of (2.1) 
then satisfies H'uII^.ooi^ir) < 4.6509 and ||-u'||j^oo(]g) < 6.4856 so that, appealing to Lemma 2.2, 
we find that any unstable spectra of the associated linearization L about u lies in the ball 
B{0,R) with 



R=l (|ir(i2)n'||i»(M) + Wf'mil^iR) +S' + 



2\2 



Next, we set Ro = 15.478 > R and ro = 1 and define Fq and Fi as in (2.14). Using 153 
evenly spaced points in A, thereby assuring relative error between successive points varies 
by less than 0.2, we compute that n(^; Fq) = for 1000 evenly spaced points in ^ in the 
interval ^ G [— 7r/6.3,7r/6.3], where here we use the scaled lifted polar coordinate method 
(described in Appendix D) to compute the Evans function D{-,^); see [BJNRZl]. 

Now, for numerical efficiency we deviate slightly from the basic numerical method de- 
scribed above and proceed to compute the Taylor expansions (2.15) of the critical modes 
near the origin. Notice that as there are generically two such branches of spectra bifurcat- 
ing from the origin at .^ = 0, there is in general no problem with making this assumption 
at this point: at the end of our numerical test, we will do a consistency check to verify 
that there were indeed only two branches to start with. Using the unsealed lifted polar 
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method (described in Appendix D) to compute the Evans function, we obtain the values 
of aj and /3j which are recorded in the q = 5 row of Table 1. Here, we have used (2.17) 
with 300 points on the A-contour, yielding relative error between successive values of aj and 
/3j less than 0.1 and 0.01, respectively, and 1000 points on the /c-contour, yielding relative 
error less than 0.2 in successive steps. From the above Taylor coefficients /?j, we can use 
Remark 2.5 to determine an appropriate low-frequency cutoff /cq- Here, we choose K = 0.5, 
R = 2 and choose /cq such that the inequality in (2.16) holds; a brief calculation shows that 
ko = 0.018365 suffices. With these choices, we use 427 points in the A-contour, assuring 
relative error between successive points varies less than 0.2 and verify that n(^; Ti) = for 
1000 evenly spaced values of ^ satisfying ko < \^\ < 7r/6.3. 

Finally we must verify that there are indeed only two spectral branches bifurcating at 
^ = from the origin. Following Step 3 of the above procedure, we must simply find 
an ri G (ro,i?o) such that n{$^;dB(0,ri)) = 2 for all |^| < k^. We note, however, that 
computing the Evans function D{\, ^) for a given ^ becomes increasingly more difficult and 
less accurate as Re(A) becomes increasingly more negative; see [BJNRZl] for more details. 
As a result, we find it more useful to compute the winding number corresponding to the 



contour 



12 



diXG 



Re(A) > -2 rnax |Re(/3j)|/cg and A - 2 max |Re(/3j)|/cg 



< 2rr 



from which we find, as expected, n(^; P) = 2 for 1000 evenly spaced values of ^ satisfying 
ICI < ^0 where again we use 427 points along the A-contour to evaluate the Evans function, 
ensuring relative error less than 0.2 between successive points. This observation, together 
with the (now justified) Taylor expansions computed above, implies that no unstable spectra 
can exist for any Bloch-frequency S,- 

From the above arguments, we conclude with great numerical certainty that the under- 
lying periodic traveling wave associated with q = 5, e = and period X = 6.3 is spectrally 
stable to small localized perturbations, in the sense of (D1)-(D3) and (H3). This procedure 
can then be repeated for different values of q, holding e = 0, the period (X=6.3) and the 
wave speed (c = 0) fixed, and the results of these computations are tabulated in Table 1. 
In particular, our numerics indicate that, as expected, there is a band of spectrally stable 
periodic traveling wave solutions, corresponding to different values of q, of (2.1) with e = 
and fixed period X = 6.3. 

Comparison with previous results. As there is a well known approximate stability 
boundary in the case e = provided in [FST], it seems worthwhile to compare the results 
of Table 1 to theirs. In [FST], the authors consider 27r-periodic solutions of the equation 



(2.211 



dtu + udxU + d^u + vd^u = 0, v > Q. 



Going to Fourier space, it is clear that when i^ > 1 all Fourier modes of the solution 
are linearly damped and hence the solution converges to a spatially homogeneous state as 



^Notice that by the previously computed Taylor expansions, we expect nim|^|<fe|, Re(Aj(5)) > 



-2maxj |7?e(/3j)|feg 



-0.0031727. 
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q 


5 


ai 


02 


/3i 


/32 


4.5 


1.41 


0.695+1.34e-05i 


-0.695+1.33e-05i 


-6.32-0.000214i 


-6.32+0.000209i 


4.6 


1.42 


6.86e-09+0.267i 


7.83e-09-0.267i 


-5.98-2.02e-07i 


-5.98-3.91e-07i 


4.7 


1.43 


-4.66e-09-0.792i 


-5.37e-09+0.792i 


-5.66-2.05e-08i 


-5.66+9.42e-08i 


4.8 


1.44 


2.44e-09-1.09i 


1.89e-09+1.09i 


-5.34-2.9e-07i 


-5.34+2.15e-07i 


4.9 


1.45 


2.38e-09-1.32i 


1.57e-09+1.32i 


-5.04-1. 15e-07i 


-5.04+8.15e-08i 


5 


1.46 


-9.44e-09-1.52i 


-9.9e-09+1.52i 


-4.74+5.02e-07i 


-4.74-2.85e-07i 


5.1 


1.47 


4.21e-09-1.7i 


3.06e-09-M.7i 


-4.45-3.9e-07i 


-4.45+2.9e-07i 


5.2 


1.48 


8.39e-10-1.86i 


2.29e-10+1.86i 


-4.17-9.99e-09i 


-4.17-7.39e-09i 


5.3 


1.49 


-1.81e-07+2i 


1.89e-07-2i 


-3.9+1.95e-05i 


-3.9+2.13e-05i 


5.4 


1.51 


-1.37e-07+2.14i 


1.67e-07-2.14i 


-3.64-4.97e-05i 


-3.64-4.72e-05i 


5.5 


1.52 


-4.29e-07-2.27i 


3.88e-07+2.27i 


-3.38-2.82e-05i 


-3.38-3.23e-05i 


5.6 


1.53 


-3.99e-09-2.4i 


-4.91e-09+2.4i 


-3.13+4.78e-07i 


-3.13-3.34e-07i 


5.7 


1.54 


2.48e-10-2.52i 


-4.01e-10-h2.52i 


-2.88-3.05e-08i 


-2.88-F2.41e-08i 


5.8 


1.55 


-7.36e-10-2.63i 


-1.31e-09+2.63i 


-2.64-hl.54e-07i 


-2.64-5.01e-07i 


5.9 


1.56 


-1.79e-09-2.74i 


-2.24e-09-F2.74i 


-2.41-h6.27e-08i 


-2.41-1. 07e-07i 


6 


1.57 


-1.08e-09-2.84i 


-1.79e-09+2.85i 


-2.18-h3.14e-09i 


-2.18+1.26e-08i 


6.1 


1.58 


3.79e-10-2.95i 


2.89e-10+2.95i 


-1.95+3.54e-07i 


-1.95-3.47e-07i 


6.2 


1.59 


1.91e-10-3.05i 


-7.21e-10+3.05i 


-1.73-4.21e-08i 


-1.73-2.22e-08i 


6.3 


1.61 


-1.14e-08-3.14i 


1.18e-08-h3.14i 


-1.52-h2.47e-07i 


-1.52-2.87e-08i 


6.4 


1.62 


-1.28e-08-3.24i 


1.37e-08+3.24i 


-1.3-h2.88e-07i 


-1.3-F2.36e-08i 


6.5 


1.63 


-3.76e-08-3.33i 


3.67e-08+3.33i 


-1.09+6.92e-07i 


-1.09-1. 71e-07i 


6.6 


1.64 


-4.03e-08-3.42i 


3.99e-08+3.42i 


-0.885+6.99e-07i 


-0.885-2.95e-07i 


6.7 


1.65 


-4.33e-08-3.51i 


4.33e-08+3.51i 


-0.681+7.4e-07i 


-0.681-2.83e-07i 


6.8 


1.66 


-2.36e-09+3.6i 


1.73e-08-3.6i 


-0.476-8.45e-05i 


-0.476-8.36e-05i 


6.9 


1.68 


-2.66e-06+3.68i 


2.66e-06-3.68i 


-0.281+0.000738i 


-0.281+0.000671i 


7 


1.69 


-1.41e-08-3.77i 


1.36e-08+3.77i 


-0.0849+1. 26e-07i 


-0.0843-2.38e-07i 


7.1 


1.7 


-1.5e-08-3.85i 


1.39e-08+3.85i 


0.109-M.6e-07i 


0.11-1.23e-07i 



Table 1: This table collects the coefficients of the Taylor expansions (2.15) of the critical 
modes bifurcating from the origin at ^ = in the case e = and X = 6.3. Here, we vary 
the integration constant q in (2.3) and note then that 6 = 6{q) is determined from q via the 
periodic boundary conditions. This table demonstrates a range of values of q corresponding 
to spectrally stable 6.3-periodic traveling wave solutions of (2.1). In these computations, 
the evaluation of the integrals (2.17) was performed using an absolute tolerance of le — 10 
and relative tolerance of le — 8 for the integration in the Floquet parameter, and an absolute 
tolerance of le — 8 with relative tolerance of le — 6 in the integration in A. Furthermore, the 
value of the Evans function along the contour was found with absolute tolerance of le — 8 
with relative tolerance of le — 10. 



t —7- oo. However, when zv < 1 there will be a finite number of excited modes leading to a non- 
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uniform solutions. Setting i^ = 1 — rj then, a stability analysis is conducted for 27r-periodic 
solutions of (2.21) with < rj < 0.69. In particular, the authors find that such solutions 
are spectrally stable for r] £ {rji^rj^) with rji G (0.25,0.3) and r/2 G (0.4,0.45) In order to 
compare these results to those provided in Table 1, we notice that by setting a = ^^"^^ — 
and introducing the new variables u = Q~^z^~^"n, x = av~'^''^x, and i = a^v~^''^ the 
27r-periodic solutions of (2.21) correspond to 6.3-periodic solutions of the equation 

dju + udxU + pdlu + d^u = 

g ^ I V . Letting v = \ — r] then, the results of [FST] suggest we should have 

stability for'p G {81,82) where 81 £ (1.32623,1.42096) and 82 G (1.65778,1.80849) in our 
yO-parametrization: these predictions are consistent with the data in Table 1 where we find 
81 ^ 1.411 and 82 ~ 1.701. 



Remark 2.6. It is worth mentioning that the numerical method of [FST] is completely 
different from ours, so that the close agreement between results gives a useful check on both 
sets of numerics. Specifically, for frequencies (^, A) bounded away from the origin, they use 
a Galerkin method similar to that used by us to generate figures using SPECTRUW; that 
is, they compute the spectra of a truncated infinite-dimensional matrix. In the numerically 
delicate small-frequency regime, they approximate the spectra by perturbation expansion 
around critical modes, as we do. However, they accomplish this in a different way, by direct 
spectral perturbation expansion of the underlying linearized operator about the wave. They 
point out also the relation between this spectral expansion and the formal modulation (i.e., 
Whitham) equations expected to govern asymptotic behavior. Though the numerics of 
[FST] are well-conditioned, and based on sound functional-analytic principles, there is no 
attempt made there to estimate numerical computation error. Indeed, for this type of 
method, this seems to us a complicated task. By contrast, Evans- function/winding number 
methods come effectively with built-in error bounds, given by the ODE tolerance and the 
requirements of Rouches Theorem to guarantee validity of the winding number. 

Second case: 7 = 1,£ = 0.2 

This choice corresponds to a generalized Kuramoto-Sivshinsky equation: the purpose here 
was to see whether we have the same qualitative picture of a band of spectrally stable 
periodic traveling waves. Here X = 6.3 is fixed. Employing our above numerical procedure, 
in an analogous fashion as for the e = case considered above, we summarize our findings 
in Table 2. As before, we find a range of values of the integration constant q associated 
with periodic traveling wave solutions of (2.1) which are spectrally stable to small localized 
perturbations, in the sense that (D1)-(D3) and (H3) hold; see also Figure 1 (b). 

2.2 A full parameter study: j = S, e^ + S^ = 1 

Finally, we conclude our spectral stability analysis by considering stability of periodic trav- 
eling wave solutions of the generalized Kuramoto-Sivashinsky equation, under the particular 
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q 


5 


ai 


02 


/3i 


/32 


4.6 


1.41 


0.694+0.2721 


-0.694+0.2721 


-6.46-1.05i 


-6.46+1.051 


4.7 


1.42 


5.44e-08+0.518i 


-6.31e-09+0.0311i 


-9.16+1.37e-06i 


-3.11-1.56e-07i 


4.8 


1.43 


-2.91e-09-0.499i 


-4.35e-09+1.05i 


-4.86+9.51e-10i 


-6.77+2.62e-07i 


4.9 


1.44 


3.51e-09+1.35i 


2.21e-09-0.791i 


-6.2+3.6e-07i 


-4.82-1.93e-07i 


5 


1.45 


-2.02e-07+1.58i 


7.16e-08-1.02i 


-5.78-1.52e-05i 


-4.64-4.23e-06i 


5.1 


1.46 


9.62e-10-1.21i 


9.17e-10+1.78i 


-4.43-6.08e-08i 


-5.42+1.06e-07i 


5.2 


1.47 


1.55e-10-1.38i 


-4.18e-10+1.96i 


-4.2+1.05e-07i 


-5.08-1.52e-07i 


5.3 


1.49 


2.81e-09-1.54i 


2.59e-09+2.12i 


-3.97-1.89e-07i 


-4.77+2.67e-07i 


5.4 


1.5 


1.43e-09-1.68i 


1.16e-09+2.27i 


-3.74-1.73e-08i 


-4.47-1.77e-08i 


5.5 


1.51 


6.84e-10-1.81i 


3.05e-10+2.41i 


-3.51-3.51e-08i 


-4.19+2.86e-08i 


5.6 


1.52 


-9.72e-ll-1.94i 


-8.41e-10+2.54i 


-3.28+2.7e-08i 


-3.92-5.11e-08i 


5.7 


1.53 


-1.05e-07+2.67i 


5.65e-08-2.06i 


-3.65+5.4e-06i 


-3.06+2. 24e-06i 


5.8 


1.54 


2.55e-09-2.17i 


2.12e-09+2.79i 


-2.84-9.26e-08i 


-3.39+6.11e-08i 


5.9 


1.55 


-2.51e-09+2.9i 


5.59e-09-2.28i 


-3.14-2.14e-06i 


-2.63-1.01e-06i 


6 


1.56 


-4.8e-10-2.38i 


-1.25e-09+3.01i 


-2.42+6. le-08i 


-2.9-5.96e-08i 


6.1 


1.57 


-4.61e-07-2.48i 


6.86e-07+3.12i 


-2.21-5.33e-05i 


-2.66-0.0001081 


6.2 


1.58 


-1.66e-09-2.58i 


1.4e-09+3.22i 


-2.01+2.63e-07i 


-2.43+1. 72e-07i 


6.3 


1.59 


-8.86e-10-2.67i 


-1.57e-09+3.32i 


-1.81+1.27e-07i 


-2.2-1. 13e-07i 


6.4 


1.6 


-4e-09+3.42i 


5.31e-10-2.76i 


-1.98+4.16e-06i 


-1.61+2.22e-06i 


6.5 


1.62 


-5.27e-10-2.85i 


-7.42e-10+3.52i 


-1.42+1.32e-07i 


-1.76+3.66e-07i 


6.6 


1.63 


1.34e-09-2.94i 


3.55e-10+3.61i 


-1.23+2.29e-07i 


-1.55-2.84e-07i 


6.7 


1.64 


-1.61e-09+3.7i 


1.98e-09-3.02i 


-1.34+5.22e-08i 


-1.05+7.35e-09i 


6.8 


1.65 


-6.87e-08+3.8i 


4.71e-08-3.1i 


-1.13+4.11e-06i 


-0. 865+1. 19e-06i 


6.9 


1.66 


-8.73e-07-3.18i 


1.39e-06+3.88i 


-0.688-3.38e-05i 


-0.925-6.86e-05i 


7 


1.67 


-9.68e-09+3.97i 


8.81e-09-3.26i 


-0.72-6.88e-06i 


-0.51-4.02e-06i 


7.1 


1.68 


-5.29e-09+4.06i 


4.25e-09-3.33i 


-0. 521+1. 23e-05i 


-0.337+6.76e-06i 


7.2 


1.7 


-2.46e-08-3.41i 


3.62e-08+4.14i 


-0.167+3.29e-07i 


-0.324-2.02e-07i 


7.3 


1.71 


-5.07e-07-3.48i 


6.83e-07+4.23i 


0.0037+0.0002031 


-0.124+0.0003591 


7.4 


1.72 


-2.72e-08-3.55i 


4.14e-08+4.31i 


0.167+1. 17e-06i 


0. 0642+1. 07e-06i 



Table 2: Similar to Table 1, this table collects the Taylor coefficients of the critical modes 
in (2.15) in the case e = 0.2 and X = 6.3. This table indicates the existence of a range of 
values of q corresponding to spectrally stable 6.3-periodic traveling wave solutions of (2.1). 
These computations were preformed with the same absolute and relative tolerances used 
for the computations discussed in Table 1. 



scaling 

(2.22) 

and e^ + (5^ = 1. It should be emphasized that the singular limit (5 — t- of (2.22) arises 
naturally in the study of small amplitude roll-waves on the surface of a viscous liquid thin 



dtu + a^(nV2) + ed^u + 5{dlu + d^u) = 0, 
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s 


Xl 


Xu 




5 


Xl 


Xu 


0.1 


9.03125 


25.7188 




0.045 


8.98047 


25.9219 


0.095 


9.03125 


25.7188 




0.04 


8.98047 


26.1406 


0.09 


9.03125 


25.7188 




0.035 


8.98047 


26.1406 


0.085 


9.03125 


25.9062 




0.03 


8.98047 


26.1406 


0.08 


9.03125 


25.9062 




0.025 


8.98047 


26.1406 


0.075 


9.03125 


25.9062 




0.02 


8.98047 


26.1406 


0.06 


8.98047 


25.9219 




0.015 


8.98047 


26.1406 


0.055 


8.98047 


25.9219 




0.01 


8.98047 


26.1406 


0.05 


8.98047 


25.9219 




0.005 


8.98047 


26.1406 



Table 3: Table of the periods Xl and Xu corresponding respectively to the lower and upper 
stability boundary of the generalized KS system, dtu + dxiu^/^) + d^u + S{d'^u + d^u) = 
as (5 ^ 0+. 



film running down an inclined plane. Indeed, one can derive (2.22) either from shallow 
water equations with friction at the bottom or free surface Navier-Stokes equations, both at 
the transition to instability of constant steady states (fluid height and velocity are constant 
in the shallow water description) for small amplitude disturbances and in the small wave 
number regime; see [W, YY] for more details on this derivation. In this case, periodic 
traveling waves solutions of (2.22) correspond to well-known hydrodynamical instabilities, 
known as roll-waves; see [NRl] and [JZN] for the spectral and nonlinear stability analysis 
of these roll-waves and for the analysis of their corresponding Whitham equations. 

In the limit (5 — )■ 0, the governing equation (2.22) reduces to the integrable KdV equation 



(2.23) 



dtu + udxU + dtu = 



where it is known that all periodic traveling wave solutions are spectrally stable to per- 
turbations in L^(M) and nonlinearly (orbitally) stable in Lpgj,([0,nX]) where X denotes 
the period of the underlying elliptic function solution; see [BD]. As discussed in detail in 
Appendix A. 3, the elliptic function solutions given in (A. 10) which satisfy the selection 
principle (A. 12) may be continued for < 5 ^ 1 to a three-parameter family of periodic 
traveling wave solutions of (2.1) parameterized by translation, period, and spatial mean 
over a period. 

Concerning the stability of these "near KdV profiles", we numerically observe that 
for (5 > sufficiently small there exist numbers Xl{5) and Xu{5) such that the associ- 
ated X-periodic traveling wave solution of (2.22) is spectrally stable to perturbations in 
L'^iM.) provided Xl{5) < X < Xu{5), with Xl « 8.49 and Xu « 26.17; see Table 3 
for details. This numerical observation complements the well known result that the "near 
KdV solitary wave" profiles of (2.22) are spectrally unstable to localized perturbations; see 
[PSU, BJNRZ2]. Analytical verification of these numerical observations, though outside 
the scope of the present work, would be a very interesting direction for future investigation. 
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See [JNRZ4] for substantial progress in this direction. 

gKS stability boundaries 




Figure 6: Plot of the stability boundaries (in the period X) versus the parameter e. Here, 
5 = Vl — e^ is fixed by the choice of e and the shaded regions correspond to spectrally 
stable periodic traveling waves. In the limits e — )■ and e — t- 1, we see the existence of 
only a single band of spectrally stable periodic traveling waves. For intermediate values, 
however, several bands emerge and the stability picture becomes much more complicated. 



At first sight, one may attempt to deduce from the previous computations 7 = 1,£ = 
0,0.2 or7 = 5— )'0,e = l that, generically, for a fixed set of parameters, there exists a 
band of spectrally stable periodic traveling waves. In fact, the situation is slightly more 
complicated and we carry out a spectral stability analysis for the full set of parameters 
e^ + (5^ = 1 and j = 6. Up to a time and space rescaling, this analysis includes all the 
particular cases treated previously: see Figure 6. This picture confirms that both for e ~ 
or (5 ~ there is a single band of spectrally stable periodic waves. However for intermediate 
values of e and 6, we clearly see that there a several bands of spectrally stable periodic 
waves and bounded bands of unstable periodic waves: this may be, among other things, 
connected to the various bifurcations of periodic waves occurring in this intermediate regime. 
Moreover, throughout this intermediate regime the nature of the transition to instability can 
be more complicated than previously seen. For example, it may happen that one diffusion 
coefficient /3j becomes positive while the other stays strictly negative; see Figure 7. This 
stands in contrast to the numerics in Tables 1 and 2. 

Comparison with previous results. The results of our numerics agree very well 
with previous results of [BaN, CDK] obtained by completely different means. In particular, 
our estimates for the limiting stability boundaries Xl and Xjj in the KdV limit are quite 
close to predictions obtained in [BaN] by formal singular perturbation analysis. Likewise, 
the global stability boundaries displayed in Figure 6 agree quite well with numerical results 
displayed in Figure 12 p. 316 of [CDK]. The methods used in [CDK] are again Galerkin- 
type as described in Remark 2.6, quite different from the Evans function methods used here. 
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Figure 7: Plot of the real part of aj and /3j from the spectral curves against the period X 
for gKS when e = 0.8 and 6 = 0.6. Figure (b) is a zoom of figure (a). 



so give a useful independent check. 

To connect our results to those of [CDK] , we give a brief lexicon between our paper and 
that of Chang et al. The figure 12 in [CDK] involves a, 5, whereas our Figure is in terms 
of parameters X and e. One has period X = 2 * vr/a and coefficient e = J/y^l + S'^)- On 
the right hand side of Figure 12 [CDK] (6 = 8/e = 0.99), corresponding to the KdV limit, 
one finds stability region a £ (0.24, 0.74) which (after translation to our coordinates) agrees 
with the bounds found in this paper and in [BaN]. The first island of instability (starting 
from e = 1) is found at 6 = 1.1 which corresponds to e = 0.74. We note that 5 = 1.1 in 
[CDK] is found to represent also the transition where new branches of periodic wave trains 
are found as S decreases below 1.1. There is also a nice discussion on the nature of the 
instabilities at this and other singularities of the stability boundary curves (p. 316 [CDK] as 
well), which we recommend to the interested reader. 



2.3 The Whitham system and time evolution studies 

In the previous section, we demonstrated numerically the existence of spectrally stable 
periodic traveling wave solutions of (2.1) in various circumstances, both in the classically 
studied case e = (Kuramoto-Sivashinsky equation) and in the case e 7^ arising from 
general thin film fiows. In each scenario considered, we found that all constant solutions 
are spectrally unstable and that, for a fixed period, an interval of the integration constant 
q (alternatively, of 5) corresponding to spectral stability; see Figure 1. 

The goal of this final subsection is to shed light on the dynamics of the underlying 
periodic traveling wave when subject to a small integrable perturbation: in particular, we 
wish to connect the observed (numerical) behavior of solutions to the theoretical predic- 
tions given by Theorem 1.1. As we will see, the long-time behavior of such solutions to 
low-frequency perturbations can be well-approximated by a formal second-order Whitham 
modulation equation, with the dynamics intimately related to the properties of these derived 
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amplitude equations. 

We consider the behavior of a fixed periodic travehng wave solution u of (2.1), corre- 
sponding, say, to a = and [3 = JS in (1.7), when subject to perturbations with characteristic 
wavenumber |i/| ^ 1, on space and time scales {X, T) = {x/v, t/v). Using a nonlinear optics 
(WKB) approximation 



(2.24) u{x, t) = uM X, T, ^^ ' ' j +uu' I X, T, ^^ ' ' 1 + 0{u 

where y >-^ u^ {X, T, y) are unknown 1-periodic functions, to find approximate solutions 
of (1.1), it follows by substitution of the ansatz (2.24) into (1.1) in the re-scaled (X, T)- 
coordinates and collecting terms of leading order that we can take 

u\X, T,y) = U{y- a{X, T) - c{P{X, T))t; (3{X, T)) , 

where a{X, 0) = and 13 {X, 0) = J3. For simplicity, we specialize now the discussion to the 
case where (3 = {k, c), with k the wave number and c the wave speed. 

As described in a more general setting in [Se], continuing the above calculation to 
higher orders, integrating over one period with respect to the fast variable y and noting 
by periodicity that integrals of perfect derivatives vanish, we find at first order in i^ the 
modulation system 

.2 25) M{k,c)t + F{k,cy, = 0, 

kt + (ck)^ = 

where k = ipx is the local frequency, c{X,T) = —ipT{X,T)/ipx{X,T) denotes the wave 
speed, M{k,c) denotes the mean of u^ over one period, and F{k,c) is the mean of f{u^). 
Here, the second equation represents simply equality of mixed partial derivatives of ■0- 

In the more specific setting of (2.1), F{k,c) is the mean of {u^)'^/2, and using the 
Galilean invariance x — t- x — ct, n — )■ u + c, we may reduce (2.25) to 

{m{k) + c)T + {H{k)-m{k)c)x=0 
kT + {ck)x = 0, 

where H{k) = F{k, 0) and m{k) denotes the mean over one period for a zero-speed wave 
of frequency k. In the classical situation, considered in [FST], where moreover e = 0, 
symmetry of the governing equation ensures m{k) = and (2.25) reduces to 

(2 27) ct + {H{k)), = 0, 

kt + (ck)^ = 0. 

Linearizing the latter about the constant (fc,0), corresponding to a background wave 
u = U{-; (k, 0)), yields the linear scalar wave equation 

(2.28) kTT - kH'{k) kxx = 
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in the local wave number k. In particular, the critical spectral curves Xj{S,) bifurcating 
from the origin at the ^ = state agree to first-order with the dispersion relation of (2.28). 
Clearly then, hyperbolicity of (2.28) is a necessary requirement for the spectral stability of 
the underlying periodic traveling wave u. 

Continuing to consider the classic case e = and a zero-speed wave, we may introduce 
similarly as in [NR2], higher order corrections to the WKB approximation (2.24) to find 
that the critical spectral curves Aj(^) actually agree to second order with the dispersion 
relation of the viscoelastic wave equation^^ 

(2.29) kTT - kH'{k) kxx = d{k) krxx 

for some d(k) depending only on the underlying wave U{.;l3). As a result, we find that 
the diffusive spectral stability conditions (D1)-(D3) are equivalent to kH'(k) > 0, cor- 
responding to hyperbolicity of the linearized first-order Whitham averaged system (2.28), 
and Re(d(A;)) < 0, corresponding to diffusivity of the second-order linearized Whitham aver- 
aged system (2.29). Thus, the periodic traveling waves below the lower stability boundary 
{5 ~ 1.411) in Table 1 correspond to a loss of hyperbolicity in the first-order Whitham 
system, while those above the upper stability boundary (6 ~ 1.701) correspond to a "back- 
ward damping" effect corresponding to the amplification of the local wave-number k{X, T) 
on time-scales of order T = t/v. These observations are discussed in more detail in [FST] 
and generalizations to cases e 7^ (and c arbitrary) may be found in [NR2] . 

Equation (2.29) recovers the formal prediction of "viscoelastic behavior" of modulated 
waves carried out in [FST] and elsewhere, or "bouncing" behavior of individual periodic cells. 
Put more concretely, (2.29) predicts that the maximum of a perturbed periodic solution 
should behave approximately like point masses connected by viscoelastic springs. However, 
we emphasize that such qualitative behavior- in particular, the fact that the modulation 
equation is of second order- does not derive only from Galilean or other invariance of the 
underlying equations, as might be suggested by early literature on the subject, but rather 
from the more general structure of conservative (i.e., divergence) form [Se, JZl].^^ Indeed, 
for any choice of /, Aj(^) may be seen to agree to second order with the dispersion relation 
for an appropriate diffusive correction of (2.25), a generalized viscoelastic wave equation. 
See [NRl, NR2] for further discussion of Whitham averaged equations and their derivation. 

We now wish to connect the formal predictions of the above WKB analysis to the 
numerically observed time-evolution of a perturbed periodic traveling wave solution of (2.1). 
Keeping with the above theme, we fix e = throughout and consider waves with period 
X = 6.3 which, after an appropriate change of coordinates, have been taken to be initially 
stationary. In Figure 8 we have fixed three periodic traveling wave solutions of (2.1). The 
wave in Figure 8(a) corresponds to a periodic traveling wave below the lower stability 
boundary so that the associated linearized Whitham averaged system formally describing 
long-time behavior fails to be locally well-posed (hyperbolic) . The resulting instability seems 



^■^The absence of terms like kxxx originates in the symmetry of the governing equation under {x, c)— > 
(-a;, -c). 

^"^ As discussed further in [Z3] , conservation of mass lies outside the usual Noetherian formulation. 
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Figure 8: Here, we present the results of four time-evolution studies for the case when e = 0. 
In Figures (a)-(c), we start with a small "Gaussian" type perturbation of the underlying 6.3- 
periodic wave (taken to be stationary by change of coordinates) and evolve the perturbation 
over time, with the vertical lines corresponding to the "peaks" and "troughs" between the 
waves. The wave-train in (a) corresponds to a wave below the lower stability boundary 
{q = 3), while the wave-train in (c) corresponds to a wave above the upper stability boundary 
{q = 10). The wave-train in (b) is spectrally stable and corresponds to q = 5.5. The "peaks" 
are plotted with thick green lines and the "troughs" are plotted with thin blue lines. In 
(a)-(d), the solid and dotted red lines indicate respectively the first order and second order 
approximations of bounds on the characteristics. They originate from a region enclosing 
three standard deviations from the mean of the perturbation. Finally, Figure (d) is the 
same as Figure (b) except for the initial perturbation is multiplied by a factor of 10. It is 
interesting to note that even when we subject the q = 5.5 wave to a large perturbation, 
lying well outside our stability theory, we observe a similar time asymptotic stability with 
analogous phase description as we developed for weak perturbations. 
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considerably more drastic than that observed for the wave in Figure 8(c), corresponding to 
a periodic travehng wave above the upper stabihty boundary. In this latter case, the first- 
order linearized Whitham system is locally well-posed but second-order diffusion coefficient, 
i.e. the coefficients /3j in (2.15), have real part with positive sign. This results in a type 
of "backward diffusion" where the amplitude of the local-wavenumber ip grows with time, 
resulting in the forced visco-elastic behavior between the individual peaks and valleys in 
Figure 8(c). For more details on the instabilities in Figure 8(c), see Figure 9. 
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Figure 9: Here, we continue the time evolution study from Figure 8(c) (corresponding to 
e = 0, X = 6.3, and q = 10). In particular, we highlight various aspects of the instability 
by zooming in on the evolution over different time intervals. Furthermore, due to the 
hyperbolicity of the first order Whitham equation (2.28) in this case the time scale on 
which the modulational instability (corresponding to Bloch frequencies |^| <^ 1) is observed 
is expected to be ~ ??~^, where t] is the difference between X and the stability boundary 
X^:] see Appendix E. 

Finally, the wave in Figure 8(b) is spectrally stable and the resulting time-evolution of 
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the perturbed wave at first sight appears to be of the form 

u{x, t) = u{x + tpix, t)) 

where ip{-,t) resembles the sum of weU-separated Gaussian waves propagating in opposite 
directions. More precisely, it appears at first sight that our small initial perturbation simply 
divides its mass into two traveling Gaussian packets which convect in opposite directions 
and satisfies limi_^oo u{x, t) = u{x) for each x G M. In particular, this description would 
suggest that 

l|n(-,t)-t2|Up(M)« ||^(-,t)||iP(R) < (l + t)-^(^-^/P\ p>2. 

However, while this intuition seems reasonable from our numerical experiments it is not 
correct: rather, the above "convecting Gaussian" description applies to the wave number ipx 
and not the phase ip. As a result, ip should roughly be described by the integral of Gaussian 
packets convecting in opposite directions, i.e. a small amplitude compactly supported (for 
each t > 0) sum of error functions of algebraically growing mass with linii^oo ip{x,t) ^ 0. 
In particular, we expect then that 

for all p > 2, and that 

It is this observation that small localized perturbations of an underlying periodic traveling 
wave solution of (2.1) behave time-asymptotically as localized shifts of the original wave 
that drives our stability analysis in the next section, even in the general case when £ ^ 0; 
see Proposition 3.2, Corollary 1, and Lemma 3.3 below. 

Finally, we remark that if one considers the more general case when e 7^ it follows 
by similar considerations that the critical spectral curves Xj{^) agree to second-order with 
the dispersion relation of a second-order Whitham modulation equation: see [NR2] for 
details of this derivation. However, the associated second order Whitham averaged system 
consists of a coupled system which, while sharing many of the properties of that of (2.26), 
is considerably more complicated to analyze directly. Nevertheless, when considering for 
definiteness again the case e = 0.2, it is found that 6.3-periodic traveling wave solutions 
below the lower stability boundary (S ~ 1.43) in Table 2 correspond to a loss of hyperbolicity 
of the first-order Whitham system (2.26) while those corresponding to waves above the upper 
stability boundary (5 ~ 1.719) correspond to a "backward damping" effect as described 
above. 

Remark 2.7. As noted previously, the modulation il^ in the above discussion is not the 
same modulation presented in Theorem 1.1, although the two are very closely related; see 
Remark 1.4. 

Remark 2.8. For rigorous justification of the Whitham equations at the nonlinear level, 
see [JNRZl]. 
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3 Proof of Theorem 1.1 

In this section, we prove Theorem 1.1, showing that spectral stabihty of a given periodic 
travehng wave solution of (1.2) implies both linearized and nonlinear stability to small 
localized perturbations. The proof closely follows the analysis of [JZl], concerning the 
analogous problem in the context of strictly parabolic second-order conservation laws. The 
main difference lies in the fact that the linear operator L defined in (1.8) is a fourth-order 
differential operator rather than second order. Nevertheless, the principles of [JZl] readily 
extend to the present higher-order case with little modification. 

3.1 Spectral preparation 

Recall that in the statement of Theorem 1.1, we normalize X to 1. Our assumptions (Hl)- 
(H2) and (D3) imply that the generalized kernel of the operator Lq, defined on L'^^j.([0, 1]), 
is of dimension 2. A crucial part of our analysis of the linearized solution operator e^* relies 
on understanding how this zero-eigenspace bifurcates from this neutral state. As expected, 
the existence of a Jordan block at A = for the operator Lq, guaranteed by hypothesis 
(H4), greatly complicates matters as compared to the degenerate case when (H4) fails; see 
[JZ3]. We thus begin our linearized stability analysis with a careful study of the Bloch 
perturbation expansion of these critical eigenvalues and associated eigen- projections near 

e = o. 

Lemma 3.1. Assuming (HI), (H2), and (D3), there exist constants ^o £ (0,7r), eo > and 
two continuous curves, j = 1,2, Xj : [—^o,S,o] -^ B{0,eQ) such that, when \^\ < S^q, 

(3.1) a{Li:)nB{0,eo) = {Xi{0,X2{0}- 

Moreover these two critical curves are differentiahle at and can he expanded as 

(3.2) A,(0= -iajC +o(0, J = 1,2 

as ^ — 7- 0. Assuming also (H3), the curves Xj are analytic in a neighborhood of ^ = 0. Thus, 
up to a possible change of ^o, there exist, for < |^| < ^O; dual right and left eigenf unctions 
{qj{6.,-)}j =1,2 o.iT'd {qj{S,,-)}j=i,2 of L^ associated with Xj{^), of form 

Qji^, •) = (iCr'l^jAO vi{^, •) + /3,, 2(6 t'2(C, •) 

qjit •) = i^PjAO M^, •) + ^2(0 V2{^, •) 
for j = 1,2, where 

• for j = 1,2, the functions Vj : [-^0,6] ^ -^pcr([0>l]) and Vj : [-Co,Co] ^ -^pcr([0, 1]) 
are analytic functions such that, when \S,\ < ^O; {''^ji^i')}j =1,2 and {vj{S^,-)}j=i^2 are 
dual bases of the total eigenspace of L^ associated with spectrum o'(L^) D B{0,eo), 
chosen to satisfy 

vi{0,-)=Ux, and i)2(0,-) = l; 
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• for j = 1,2 and k = 1,2, the functions f3j^k '■ [— ■^o^'^o] ~^ C and /3j^k ■ [— Coi'^o] ~^ C 
are analytic. 

Finally, assuming in addition that (D2) holds, the spectral curves Xj can be expanded as 
(3.3) x,{0 = -iaj^-b,e + Oi\C\''), j = l,2 

for some aj G M and bj > 0. 

Proof, (following [JZl, JZN, NR2]) First, as already mentioned in Remark 1.3, assumptions 
(HI), (H2) and (D3) ensures the possibility of a parametrization /? = (k, M) by wave number 
and mean We adopt such a parametrization. This provides LqUx = and LodcU^-; f3) = 
-dMc{i3)ux with {v,Ux)lI^^([o,i]) = and (i), 9m^(-;^))l2_^^{[o,i]) = 1> where v = I. 

Since L^ has dense domain H^^^{[0, 1]) compactly embedded in Lpgj.([0, 1]), its spectrum 
consists of isolated eigenvalues of finite multiplicity [He] . As is separated from the rest 
of the spectrum of Lq, assumption (D3) and standard spectral theory for perturbations by 
relatively compact operators (see [K]) yields constants ^o £ (0, vr), eq > and continuous 
functions Ai, A2 such that, for |^| < ^0, ^(-^g) ni3(0,eo) = {^i(?)) ^^2(0}- Moreover this 
also yields analytic dual right and left spectral projectors associated to spectrum in B{0, eq). 
Analytic dual bases of the right and left eigenspaces may then be obtained by projecting 
dual bases for spectral spaces of the spectrum of Lq in B{0,eo). By the first paragraph of 
the proof and the conservation law structure of the governing equation, we may choose such 
bases in the form {ux, dMU{-; /?)} and {*, v} and obtain in this way the functions {vj} and 
{vj} of the lemma. 

We have now reduced the infinite dimensional spectral perturbation problem for the 
operators L^ to the spectral analysis of 



^i = [{'^ji^r),L^vi{t-))Lmo,i])_ 



j,i 



a 2 X 2 matrix perturbation problem. By direct calculation Mq = [ ^ n ) ' Below, 

however, we will scale M^ to blow up at the double eigenvalue. To do so, we expand the 
operator L^ as 

(3.4) h = Lo + ^e^(i) + m'L^^^ + {iCfL^^^ + m'L(^^ 
and note, specifically, that 

(3.5) lW = -(n-c) -3ea2- 2(55^ -49^. 

By either a direct calculation or by first scaling the parameterization then differentiating the 
profile equation with respect to k (see [NRl, NR2]), we find that ('i'2(0, ■), L^^'vi{0, •)) = 0. 
Therefore, M^ can be expanded as ^ — )■ as 
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so that, in particular, the scaHng 



(3.6) M^ := iiO-'S{OM^S{0-\ S{0 := ( 



i^ 
1 



preserves smoothness in ^ at ^ = 0. Since the eigenvalues Tnj{^) of M^ are (^0~"'^^i(0) their 
continuity implies the differentiability of the functions Aj(^) at ^ = 0. Assuming, in addition 
to above, that assumption (H3) holds, it follows from the fact that mj{0) = —aj that the 
eigenmodes of M^ are analytic in ^ in a neighborhood of ^ = 0. Undoing the scaling finishes 
the proof, up to the observation that (3.3) follows from (D2) and Xj{S,) = Xj{—S,). D 

3.2 Linearized bounds 

We begin our stability analysis by deriving decay rates for the linearized solution operator 
e of the linearized equation (1.8). Recalling the inverse Bloch transform representation 
(1.13) of the linearized solution operator, we first introduce, as in the proof of Lemma 3.1, 
for each S, G (— Co,Co) the spectral projection P{^), analytic in ^, onto the total eigenspace 
associated with spectrum cr{L^) n 5(0, eo) = {'^i(O) ^^2(0} of the Bloch operator L^. We 
also choose a smooth cutoff function (/) : [— 7r,7r) — )■ [0, 1] such that 

,,,. _ r 1, if lei < ^o/2 

and split the solution operator S{t) := e into its low Floquet-number critical part 

(3.7) S'{t)gix) := f e*«"0(O[P(Oe'^«*5(e, ■)Kx)dC 
and exponentially-stable part 

(3.8) S''it)gix) := I e*«-[(l - H^PiO) e^^'gi^, .)]ix)dt 



— IT 



As the analysis of the critical part is considerably more delicate, we being by deriving L^ 
bounds on S^^{t). By standard sectorial bounds [Pa, He], the fact that L^ and H^ spectra 
coincide and the spectral separation of the Xj{^) from the remaining spectrum of L^ we 
have bounds 

(3.9) ||i^re^«*(l- 0(6^(0) 5llH,^„([o,i]),«<Ct-"^e-^*||<7l|^;,,([o,i]),5 

for some constants 9,C > 0, with ||5'll^r^^([o,i]),^ = Z]j=o Wi^^ + ^0^ dWi^ ([o,i])- Using the 
fact that L^ is a relatively compact perturbation of the fourth order differential operator 
(—1 — {dx + ^C)^)) i^ conjunction with (1.14), we immediately have the H^ bounds of the 
following Proposition. 
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Proposition 3.1. Under assumptions (H1)-(H4) and (D1)-(D3), there exist constants 
C,9 > such that for all 2 < p < oo, < Ali + I2 = k < K + I, < Ami + m2 + m^<K, 
K as in (HI), r = 0, 1 and t > 



d['d'iS'\t)d'ig 



. 4ii+i2+i3 _ 

< Ct 4 e Wr 



lo^ a,j. o \T-)o^ y rpcin,-) ^ <-^i ^ ^' e ||y||L2( 



1 (\ \\ 4m]^+m2+m3 

< r.-t 

\LV{ 

\S'\t)g\l,,^,<Ce-'%\\HH 



\LP( 

Proof. By the above discussion, all that is left is to explain how to verify the stated L^ 
bounds. These follow from L^ and H^ bounds and the Sobolev embedding inequality 

Next, we analyze the critical part of the solution operator S{t). For this purpose, it is 
convenient to introduce the (critical) Green kernel 

G\x,t-y) := S\t)5y{x) 

associated with S^ and 

[Gj(.,t;y)](x):=0(e)P(Oe^«*NW 

the corresponding integral kernel appearing within the Bloch- Fourier representation of G^, 
where [•] denotes the 1-periodic extension of the given function onto the whole real line. We 
first prove the following lemma yielding a spectral representation of G^ incorporating the 
results of Lemma 3.1. 

Lemma 3.2. Under the hypothesis (H1)-(H4) and (D1)-(D3), we have 

2 



[G\{-,t;y)] {x) = 0(6 j;e^^«\-(^,x)g,(e,y)*, 



i=i 



(3.10) G\x, t;y)= e*«(--^) [gJ(-, t; y)] {x)d^ 



2 

e^^^--yU{OY.^^'^^^''^jitx)qj{tyrd^, 
i=i 

where * denotes the matrix adjoint, or complex conjugate transpose, and qj{^, •) and qj{$., ■) 
are right and left eigenf unctions of L^ associated with the eigenvalues Xj{C) defined in (3.1), 
normalized so that {qj{^, •), qj{^, ')) l'^ (\o ii) = 1- 

Proof. The first equality follows by the spectral decomposition of e 5*^ and the spectral 
description of Lemma 3.1. The second equality follows by the inverse Bloch transform 
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formula (3.7) and Fourier transform manipulations using the fact that both discrete and 
continuous transforms of the centered (5- function are constant equal to {2tt)^^ to get 



Syitx) 



1 

2^ 



lei 



2inl{x—y) 



-i^y 



[6y]{x). 



The third quality now follows by substitution; see [0Z4] for further details. 



D 



Continuing, we point out that it seems unlikely that the low-frequency Green function 
G will satisfy L^ — )■ L^ bounds which are suitable for our purposes. To see this, notice 
from Lemma 3.2 and Lemma 3.1 that we have the representation 



G'{x,t;y)= / e*«(-^^)0(e)j;e^:'«)* (ie)-Xi(m,2(eK(e,:r)^2(e,?/r + 0(l) d^ 



of the critical Green kernel. From assumption (D2) we expect, for example. 



\G'{;t;y)\ 



L°=( 



C^\C\-'e''^^^\{0 



Li([-^,^)) 



where the right hand side is interpreted in the principal value sense, which is merely bounded 
and hence does not decay in time. In order to compensate for this lack of decay arising 
from the Jordan block associated with the translation mode at ^ = 0, we separate out the 
bounded translation mode from the faster-decaying "good" part of the Green kernel. To 
this end, notice that by defining the function 



(3.11) 
we have 



(3.12) 






G^{x,t;y) - u^{x)e{x,t;y) 



M-y)^(^\ V^ .A,«)t/,. ^ ^,^ ^i(e,x)-7;i(0,x) . 



+ 



V(oE^- 



%,i(e)- 



QjiCyTdC 



Now since all quantities involved in (3.12) are C in ^, from (D2) we deduce 



(3.13) \\G'i-,t;y)-u^i-)ei-,t;y)\ 



<C 



^ ^ e 



'0\(\H 



HO 



LH[-n,n)) 



<G{i + ty 



'1/2 



yielding algebraic decay of G — UxS in time. The more general effect of this regularization 
is the content of the following proposition. 
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Proposition 3.2. Under the assumptions (H1)-(H4) and (D1)-(D3), the critical Green 
kernel G^{x,t;y) of (1.8) may be decomposed as 

G\x, t; y) = u^{x)e{x, t; y) + G\x, t; y) 

where for all t > 0, l<q<2<p< co, and 1 < r < 4 the residual G {x, t; y) satisfies 

&{;t;y)g{y)dy 



<C{l + t) H-? py 



\9\\li( 



(3.14) 



d;G'{;t;y)g{y)dy 
dt&{;t;y)g{y)dy 



LP( 



LP( 



LP( 



<C{l + t) Hi "p) 5||^||^ 



<C{l + t) Hi p) '^WgU.^^y 



Furthermore, for all t >0, l<q<2<p< oo, < j, /, j + / < K, and 1 < r < 4 we have 



(3.15) 



didldle{-,t;y)g{y)dy 



<c{i + ty 



ifi i\ u+n 



2\q V 



mm 



LP( 



and, if in addition j + I > 1, or (j = I = 0, p = oo and q = 1), 



(3.16) 



didte{-,t;y)g{y)dy 



<C(l+t)2 2V9 pj 2 



9\\li( 



LP( 



Proof. To begin, let G^ {x, t; y) := G^ {x, t; y) — Ux{x)e{x, t; y) where the function e is defined 
in (3.11). By interpolation, it is sufficient to consider only the cases ^ = 1 and q = 2. 
(i) Case q = 1. To prove bounds for g = 1 we rely on Triangle Inequality 



F{-,t;y)g{y)dy 



LP( 



< ll5'llLi(R)Sup||F(-,t;y)||ip( 



Now, to generalize (3.13) and prove the first part of (3.14) in the case q = 1, we use (1.15) 
to get for any {t, y) 



&{;t;y) 



<C 



LP( 



^ ^ e-''i«i'V(e) 



LP' ([-TT,Tr)) 



<C{l + t) 



|(i-i/p) 



with p' the Holder conjugate of p, 1/p + 1/p' = 1. The second part of (3.14) comes in a 
similar way by observing that for j = 1, 2 and r > 1, the fact that Qj(0, y) = /3j,2(0) implies 



(3.17) 
and thus 

d;&{;t;y) 



dlqM,y) = d;q,{^,y)-d;q,{0,y) = O{0, 



< C 



LP( 



c^ieie-'|«|'V(o 



LP' ([-TT,Tr)) 



<c{i + t)-^2^'-'M-l 
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Likewise the last part of (3.14) stems from the fact that the time derivative brings down 
factors XjiO = 0(0- 

Using again (3.17), we obtain for r > 1 



did^,d;e{x,t;y)= / e*«("-^)</.(0 J] eV«)*(iO^A,(0'/?/,i(0 

i'=i 



d;qr{^,yr-d;q,,{o,yy 



< 



and another use of (1.15) provides 



dididle{;t;y) 



LP( 



<c 



e ^ ler'e-^i^i v(6 



LP' ([-7r,7r)) 



<c(i + t)- 



■i{i-i/p)-|(j+0 



This proves (3.15) in the case g = 1. Inequality (3.16) follows similarly when j + / > 1. 

Finally, we must obtain (3.16) in the critical case j = / = 0, p = oo and g = 1. Recalling 
(3.3) we introduce A(^) := —iaj^ — bj^"^ and expand 



e{x,t;y)= / e*«(""^),/.(0 |^ e^ 



xAOt^jAOQjity)* - ^j,mqjio,yr 



i=i 



iC 



di 



+ 



2 



Mo 



^e(^:'(«)-^^(«))*-l 



(3.18) 



+ / e 

'-TT 



iii{x-y)^^\M)t <t^^^)-^ 



i=l 









/3,,i(0)42(0)d^ 



e*«(-y) ^ e^^«)*^/3,- i(0)42(0)de 



2 



.r 



1 



+ P.V. / e^«("-^)ye^:'(«)*^/3,-i(0)/3j-2(0)de 

All but the last term in (3.18) may be bounded using Haussdorff- Young estimates, either 
using the classical one or (1.15). To complete the proof then, we must derive an appropriate 
L°^ bound on the last integral in (3.18). Since 



1 



Mix-y)^HOt 



di 



^-{x-y-ajt)^/4bjt 



'2'^ Jr" " " y/47rbjt 

the final principal value integral in (3.18) is recognized to be 

2vr^/3,,i(0)/3,,2(0)errfn('^^|=^ 



J = 1,2, 



hence is bounded in L°^ as claimed. 
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(ii) Case q = 2. To prove bounds with q = 2 we directly apply (1.15) to Bloch formu- 
lations. For instance expand 



&{x,t;y)g{y)dy 






J^x 



0(OEe^^«)*/3,,i(O 



i=i 



t;i(g,x) -fi(0,a;) 



{Qji^, ■), gi^, ■)) Ll^J[o,i])d^- 



The generalized HaussdorfF- Young estimate (1.15) thus yields 



G\;t;y)giy)dy 



LP( 



< c 



<c 






< C(l + t) H2 P 



LP ([-tt.tt]) 
S'llLafR) 



where p' is such that l/p+ 1/p' = 1, and p" such that l/p' = 1/2 + 1/p". The other bounds 
on & follows in the same way, observing that space derivatives dy allow the use of (3.17) 
and thus provide an extra 0{^) factor in the appropriate integrals. The bounds on e can 
be obtained in an analogous way, completing the proof. D 

Finally, we combine the above various exponentially-stable and critical bounds to obtain 
decay estimates on the Green function 

G{x,t;y) = S{t)6y{x) 

associated with the full solution operator S{t) = e^^ . Prior to that, we let x ■ [0) oo) — t- [0, 1] 
be a smooth real valued cutoff function such that 



xit) 



0, if < t < 1 

1, if t > 2 



and define 
(3.19) 



e{x,t;y) := xit)eix,t;y); 



the purpose of the time cutoff function x will be made clear in Remark 3.4 below. The 
following result immediately follows from Lemma 3.2 and Proposition 3.2. 

Corollary 1. Under assumptions (Hl)-(H^) and (D1)-(D3), the Green function G{x,t;y) 
of (1.8) decomposes as 



Gix, t; y) = na;(x)e(x, t; y) + G{x, t; y) 
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where, for some C,9 > and all t > 0, 2 < p < oo and 1 < r < 4, 
(3.20) 



d:^G{;t;y)g{y)dy 



< mill 



(3.21) 



G{-,t;y)g{y)dy 



LP( 



LP( 



Ct 4V 2 pj 4 



+ C(l+t) H2 p) 5 



\9\\lH 



-if 3_l\ 1 I r 

(1 + t) H2 Py* 2 + 4 



IS'IIli 



< C(l+t) H^ P. 



l5'llLi(IR)nHi(]R) 



Furthermore, e{x, t;y) = for < t < 1 and there exists a constant C > such that for all 
t>0, l<q<2<p<oo,0<j,l,j + l<K, and 1 < r < 4 we have 



(3.22) 



didtdle{-,t;y)g{y)dy 



<C{l+t) H9 Py 



(j+fc) 
2 



\9\\li( 



LP( 



and, if in addition j + I > 1, or (j = I = 0, p = 00 and q = 1), 



(3.23) 



didle{-,t;y)g{y)dy 



<C{l + t) 



1 1/1 i\ (j+fc) 



2 2 \(j p 



IffllMf 



LP( 



Remark 3.3. As may be clear from the proofs of various linear estimates discussed above, 
except for the critical bounds on e or e in the case j = l = 0, p = oo and q = 1, the in- 
troduction of the critical Green kernel was a pure presentation device. Indeed, we adopted 
here it to mark proximity with [JZl]. Yet, since almost all the above linear bounds were 
proved using Haussdorff- Young type estimates, there may be some gain in clarity and ef- 
ficiency in presentation in keeping all descriptions at the level of semigroups and Bloch 
symbols. This latter approach was recently adopted in [JNRZl, JNRZ2, JNRZ3], where 
further decompositions of the critical part of the solution operator are needed. 

3.3 Nonlinear preparations 

Given the linearized bounds on the linear solution operator S{t) = e derived in the 
previous section, we are now in position to consider the effect of the small nonlinear terms 
that were omitted in obtaining the linearized equation (1.8). Our first task is to explain 
how to implement at the nonlinear level the separation of critical phase-shift contribution. 
To this end, let u{x, t) be a solution of 



(3.24) 



dtu — cdxU + d^u + ed^u + Sd^u + dxf{u) = 



and define the spatially modulated function n(x, t) := u{x+ij){x, t), t), where ip : M.xM.^ — t- M 
is a function to be determined later. Moreover, let u{x) be a stationary periodic solution of 
(3.24) and define the nonlinear perturbation function 



(3.25) 



v{x, t) := u{x, t) — u{x). 
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Lemma 3.3. For v as above, we have 

(3.26) i^t-L){v-^Pu,)=M, N = d^Q + d^U + dtS 

where 

Q = - {f{u + v)- f{u) - df{u)v) 

n = 1ptV - d——^Vx - £-— —Ox -— —Vx - £Ox 



l + 1px l + 1px \l + lpx J \l + 1px 

^ -9x ( , . , Ox ( , , Vx ]] - dxi ——^dx ( , . ^ Vx 



l + lpx \l + 'lpx \l + 1px J J \l + 1px \l + i). 

pa. ( -i'x \ ^ i^l _ -i^x p. ( -iJx -\ i^l - 

dx rr—Vx - b———Ux - e———Ox , . , u^ - g , , , u.. 



l+1px J l + lpx 1 + V'x Vl + V'a; / 1+V': 



XX 



Jx 



--^^ [St-a - T^^^ (t^^^ {^A^ - ^^i (^-^ 

\i + iJx J i + i^x \i + il^x \i + ipx J J 1 + il^x \i + i^x 

~Vxxx C>x \ ~ \ ~j~Ox \ ~ \ ~J~Ux ] ] Ox [ ~ ; ~J~Uxx I Ox 



1 + ^Px ^\i + i^x \i + i^x JJ \i + i^x J "Vi + ^^ 

S = -vil^x- 
Proof. From the definition of u above, it follows that 

(1 + iJx)dtu - {11)1 + c)dxu + dxf{u) + 5dx ( ———dxU 
+ edx ( —dx ( —dxU 

Hence, subtracting the profile equation —cux + Uxxxx + ^Uxxx + ^Uxx + {f{u))x = implies 
{dt - L)v - {'ipux)t = dxQ- ipxVt + iptVx - 5dx ——^{u + v)^ 

\l+'lpx 

Vi + Vx yi + V'x J J ■" \i + 4^x 

— dx I ^dx ( —dx ( —(u + v)x 

.2f-i'Xr,f 1 .., , . \\ P,Z( -V'X 



''Vi + V'x Vi + V'x yy "^Vi + V-x 

or, equivalently, 

(9i - L)v - {ipux)t = dxQ + dxTZ + dtS - 5dxi-ipxUx) - edx{-^xUxx) - edl{-ipxUx) 

- dxi-'lpxUxxx) - dl{-1pxUxx) - dl{-1pxUx). 

Using the fact that tpLux = completes the proof. D 
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Using Lemma 3.3 and applying Duhamel's principle, recalling that e{x, t;y) = for 
< t < 1, we obtain the implicit integral representation 

v{x,t) = Ux{x)ij{x,t)+ / G{x,t;y)v{y,0)dy+ / / G{x,t - s;y)Af{y,s)dy ds 

Jr Jo Jr 

for the nonlinear perturbation variable v and some function ijj still to be determined. Re- 
calling Corollary 1, it follows that by defining ip implicitly via the integral formula 

(3.27) '4'{x,t) = - / e{x,t;y)v{y,0)dy- / / e{x,t - s;y)Miy,s)dyds 

Jr Jo Jr 

the perturbation variable v must satisfy the integral equation 

(3.28) vix, t)= [ Gix, t; y)v{y, 0)dy +11 G{x, t - s; y)M{y, s)dy ds. 

Jr Jo Jr 

Furthermore, recalling that e{x,s;y) = for < s < 1, we find by differentiating (3.27) 
that 

(3.29) did^^tPix, t) = - f did^eix, t; y)v{y, 0)dy -If d{d^e{x, t - s; y)M{y, s)dy ds. 

Jr Jo Jr 

for < i < 1 and < /c < K + 1. 

Remark 3.4. The purpose of the time cutoff function x introduced in (3.19) is precisely 
to ensure that ip{-,0) = 0, corresponding to the fact that we are dealing with perturbations 
that are initially localized in space. Yet, in this periodic context it is reasonable to consider 
also perturbations which affect the phase or wave number of the underlying periodic profile, 
which corresponds to perturbations which are not initially localized in space. Stability of 
periodic wave trains to such non-localized perturbations have been the focus of much recent 
study; see [SSSU, JNRZ2, JNRZ3] for analysis in the context of reaction diffusion systems, 
and [JNRZl] for analysis in systems with a conservative structure. 

Remark 3.5. As discussed in Remark 1.4, the implicit ^-dependent change of variables 
performed above has the effect of introducing at the linear level the critical non-decaying 
translational mode ipUx, while at the same time ensuring that only derivatives of ip, which 
decay in L^(M), appear in the nonlinear terms A^ in Lemma 3.3, as opposed to non-decaying 
terms involving ip itself. 

Equations (3.28) and (3.29) together form a complete system in the variables (f , V'tjV'a;) 
and from the solution of this system, if it exists, we may recover the phase modulation 
tp through (3.27). Furthermore, the short-time existence and continuity with respect to t 
of solutions {v,tlJt,^Px) e H^{R) x H'^{R) x H^+'^{R) follows from (3.26) together with 
(3.29) and a standard contraction-mapping argument based on (3.27), (3.22), (3.23), and 
the following nonlinear damping estimate. 
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Proposition 3.4. Assuming (HI), there exist positive constants 6, C and eq such that if v 
and ip solve (3.26) on [0, T] for some T > and 

sup \\{v,^x){t)\\HKm+ sup ||Vt(t)||j:/if-i(R) < eo 
te[o,r] te[o,T] 



then, for all < t <T, 



k(t)||^K«, <Ce-^*||<0)||2 



(3.30) /•* , ^ / \ 

+ C / e-^(*-^) (lb(5)lli2(j,) + ||Vx(5)||^K+l(M) + ||Vt(5)ll?,K-2(M)) ^5. 

Proof. First rewrite (3.26) as 
■Ui + S^-u = -e9^v-(55^i; + c9:Ei;-9a;(/(w + f) - /(ii)) - -— -^5:E(/(n + i;)) 

1 + V'x 

+ -I , , — [Ux + Vx) -0 dx , , . {ux + Vx) ] -odxl —-—-{ux + Vx, 

(^x ( , . , dx ( , . , {ux + Vx)]] - edx ( ——^dx ( , , , {ux + -v^^ 



1 + Vx Vl + V'x Vl + V':. JJ ^l + Vx Vl + V': 

— ea^ ( ——^{ux + Vx)] - ——^dx ( ^ , , dx ( , , , ^x- ( , , , (^x + Vx] 

Vi + V-x 7 i + V'x \i + i^x \i + i^x VI + ^^0; 

— edx I ^dx I —dx I — (^^x + Vx] 

— e9^ I —dx ( —{ux + Vx) I I — e5^ ( 1— (^x + Vx) ] . 

By taking scalar product against Y^j=Q{—^ydiv and using integration by parts and Sobolev 
embedding Hfi'llLoom-) < C||5f||//W]g), we obtain under a smallness assumption as in the state- 
ment of the proposition 



I d / 2 ^ ' ^ 1 



) it) + ^\\vmlK+2^^) < C(||i;(t)||^K+l(K) + ll^x(t)|lWl(R) + ll^t(0llHi^-2(M)) 

Using now the Sobolev inequality 

with a sufficiently small rj reduces the result to a simple integration. D 

3.4 Nonlinear iteration 

With the above preparations in hand, we are now prepared to state the main technical 
lemma leading to the proof of Theorem 1.1. For this purpose, associated with the solution 
{v,iljt,tpx) of the integral system (3.28) and (3.29) considered in the previous section, define 

(3.31) v{t) ■= sup \\{v,iJt,'4^x,ipxx){s)\\HKm{l + s)^/*. 

0<s<t 
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By standard short-time i7 (R) existence theory, 7] is continuous so long as it remains 
sufficiently small. Using the linearized estimates of Section 3.2 we now prove that if r/(0) is 
small then r]{t) remains small for all t > 0. 

Lemma 3.6. Under assumptions (Hl)-(H^) and (D1)-(D3), there exist positive constants 
C and £0 such that if v{0) is such that 

Eo ■■= \\v{0)\\Li(R)nH'<{K) < eo and r]{T) < eo 
for some T > 0, then for all < t < T we have 

(3.32) 7]{t)<C{Eo + r]{tf). 

Proof. First, note that by Lemma 3.3, we have under the smallness assumption on t] that 

(3.33) \\iQ,T^,Smh^u)nLHR) < C7?(t)2(l + t)-V2 

for some constant C > 0. The contribution of dfS needs some special care. For this 
purpose, note that, integrating by parts in critical contribution, we obtain from previous 
linear bounds and 5(0) = 

ft 



f f G{;t-s;y)dsS{y,s)dyds < C||5(t)|U2(M) + C / e'''^'-'^\\dtS{s)\\HiiR)ds 

Jo Jm. lp{r) jo 

[ (l+t-s)"H^"p)"5||5(g 

Jo 



+ c I a + t-s) n^ -pj -^\\s(s)\\ru^,ds. 



Now, observe that, using (3.26) to bound vt, we have under the same smallness assumption 
the estimate 



(3.34) \\SmHum<Cvit)H'^ + t) 



-1/2 



Using then the second part of (3.20) and (3.21) to bound other terms, we obtain from (3.28) 
the estimate 

\\vm^,^^^<Cil + t)-'^i'-p){Eo + r^{tf) 

(3.35) +Cri{tf (t-s) ns p) -*(i + t-s) ^5 p) ^{l + s)-^ds 

Jo 

<Cp{Eo + i^{tf){l + t)-'^i'-p) 
and similarly, using (3.22) and (3.23) with g = 1, we obtain from (3.29) the estimate 

ll(^t,V'x)llvF^+i.P(R) < cii + ty'^'^^'-p^Eo 

(3.36) + Cr,{tf [ {l + t- s)-^(i-i/p)-| (1 + s)-i/2rfs 

Jo 

< Cp {Eo + r,{tf) (1 + ty 



PJ 
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valid for all^^ 2 < p < oo. The assumed smallness guarantees that we may apply Proposi- 
tion 3.4 and get 

lb(t)||^K(jj) < Ce'^'Ei + C{Eo + vitff J^ e-^(*-^)(l + sy'^'ds 

Together with (3.36), this completes the proof. D 

Proof of Theorem 1.1. We are free to assume the constant C in (3.32) is larger than 1. 
Since rj is continuous and r/(0) = ||t'(0)||j|/K(]g) < Eq, it follows by continuous induction that 
if AC'^Eo < 1 then r?(t) < 2CEo for ah t > 0. We may then use (3.35) and (3.36) to get 
uniform bounds for p G [2,4] 



lb(i)llLp{R) + ll(^t,V'x)(i)lliyK+i,p(K)<C^o(l + i) ^^' ^ 
Next, we rewrite TZ in form 

7^ = 7^l + a^7^2 + dln3 + d^n^, 

+ Tlf'{x, '4)x{x, t),1pxx{x, t), i^xxxix, t),1pxxxxix, t))v{x, t) 

so that 

||(Q,7^l,7^2,7^3,7^4,5)(t)||i2(K) < C(||7;(t)||i4(K) + ||(V't,V'.)(t)fv^K+M(ij)) < Ci?o(i + t)-^ 
and 

\\{dxQ,dxni,dln2,dln3,d^n4,dtS){t)\\^,^^^ < CEo{i + t)-h 

Using now the estimate 

ft 



Jo 

+ C /(l + f-s)"K^-?)-^||5(5)||i2(M)d5, 



t 

G{-,t - s; y)dsS{y, s)dy ds 





together with the first part of (3.20) and (3.21), we obtain for all p G [2, oo] 

\Ht)\\LviM.)<C{l + t)~'^^^^^)EQ + CE(, f (l + t-s)"K^-^)-5(i + s)-i(i5 

Jo 



< CEo{l + t) n^ p 



^^Notice that the above integral estimates fail in the case p = cx3 due to a term of size log(l + t) arising 
from integrating on [|, i] • 
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and, similarly, using (3.22) with q = 2 and (3.23) both with q = I and q = 2 we obtain the 
estimate 

IKV't, V'x)(t)||vK^+l,P(R) < C(l + t)"Kl-p)^o + CEo f {1 + t- 5)-K5-^)-5 (1 + ^)-| d, 

Jo 

This completes the proof of (1.20). 

To establish (1.21), note that (3.22) and (3.23) imply 

\mt)\\L^iR) < CEo + CEi J (1 + t - s)"K^-p) (1 + s)-'2ds < CEo. 

Finally, notice that by definition we have that 

u{x, t) — u{x) = v{x, t) + {u{x, t) — u{x + ip{x, t), t)) , 

hence 

\\u{t) - u||loo(i;) < ||t;(t)||i^(]R) + ||na;||Loo([o,i])||^(t)||Loo(K) < CEq. 

We thus obtain the final L^ H H^ — t- L°^ bound (1.21) of Theorem 1.1. This completes 
the proof of the main theorem, hence establishing the nonlinear L^ n H — )■ L°^ stability 
of the underlying periodic traveling wave u under the structural and spectral assumptions 
(H1)-(H4) and (D1)-(D3). 

Ackno'wledgement: The numerical stability computations in this paper were carried 
out using the STABLAB package developed by Jeffrey Humpherys and the first and last 
authors and the SpectrUW package developed by Bernard Deconinck and collaborators; 
see [BHZ2, CDKK] for documentation. We gratefully acknowledge their contribution. We 
thank also Indiana University Information Technology Service for the use of the Quarry 
computer with which some of our computations were carried out. 

A Appendix: Survey of existence theory for generalized KS 



In this appendix we give a brief survey of the existence theory for periodic traveling wave 
solutions of the generalized Kuramoto-Sivashinsky equation (1.1) with nonlinearity f{u) = 
3u^. For the forthcoming discussion we will also consider the particular scaling where 7 = 5, 
i.e. the equation 

(A.l) ut + 6uux + ed^u + S{dlu + d^u) = 0. 

Under this scaling, there are three situations to consider: the rather generic case where 
£ = 1 and 6 = 0(1); the classic Kuramoto-Sivashinsky limit 5 = 1 and |e| <^ 1 and the 
(integrable) Korteweg-de Vries (KdV) limit e = 1 and < 5 <^ 1. In the first two cases, 
we obtain only existence of small amplitude periodic wave trains through a normal form 
approach whereas we obtain existence of large amplitude periodic traveling waves in the 
KdV limit, by using Fenichel theory. 



A APPENDIX: SURVEY OF EXISTENCE THEORY FOR GENERALIZED KS 54 



A.l The generic case: Hopf bifurcation analysis 

Traveling waves u{x,t) = U{x — ct) of (A.l) are then readily seen to be solutions of the 
ODE 

(A.2) 6{U"' + U') + eU" + 3U^ - cU = q, 

where g G M is a constant of integration. We begin by considering the generic situation 
where e > and 5 j^ are fixed and arbitrary. In this case, so long as c^ + 12q > 0, (A.2) 
possess two stationary solutions U-{c,q) < U+{c,q) such that 3U^ — cU± = q. Linearizing 
(A.2) about a constant state U = U± thus yields 

(A.3) (6C/± - c)U + eU" + 6{U' + U'") = 0, 

which is seen to undergo a Hopf bifurcation (necessarily) at U+ when CHopf (?) = 6C/+ — e. 
Denoting V = U — [/+ and 6C/+ — c = e + fi with |^| <^ 1, equation (A.2) then reads 

(A.4) {3V + e + fi)V + eV" + SiV' + V'") = 0, 

from which we find, for sufficiently small ^ on one side of zero, a family V{iJ,) of periodic 
orbits with amplitude Od/up^) and frequency k{fi), with /c(0) = 1 and fc'(O) 7^ 0. As a 
result, in the neighborhood of the Hopf bifurcation we find a three dimensional manifold of 
small amplitude periodic traveling waves parameterized by translation, the wavenumber k, 
and the integration constant q. 

Remark A.l. Note that when c^ + 12q = 0, the equilibrium states coincide. The nature 
of the bifurcation occurring at this degenerate point requires a more delicate normal form 
analysis, which we postpone to the next section. 

Finally, notice that for |/i| <^ 1 the solution U of (A.2) satisfies 

where {U) denotes the spatial mean of U over a period, so that, by the definition of q, we 
find 

3{Uf-c{U) = q + OM'/^). 

Setting M{k, q) = (U) it follows that dqM is non-zero by the definition of /U. Here, we are 
using that e > is fixed. In particular, in a neighborhood of the Hopf bifurcation, we can 
switch from the {k, q) parametrization of the local manifold of periodic traveling waves (here 
identified up to translation) to the {k,M) parametrization; see Remark 1.3. This provides 
the following proposition. 

Proposition A.2. Let M gM. be fixed and k < 1 be such that 1 — k is small. Then there 
exist a unique c{k,M) and a unique q{k,M) such that there exists a 1-periodic U{- 'jkjM) 
solution of 

(A.5) k (3[/2 - cU) + ek'^U" + 6 {k'^U' + k^U'") = q, (U) = M . 

Moreover this solution is unique up to translation. 
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A. 2 The classic KS limit: a normal form analysis 

Here, we complement the above bifurcation analysis in the previous section with a normal 
form analysis near the degenerate case |e| <^ 1. In order to find periodic solutions, we 
necessarily work in the neighborhood of < c^ + 12q <^ 1. Setting V = U — [/+ then we see 
that V must satisfy (A. 4) which we can rewrite here as the first order differential system 
(A.6) 

x' \ / 1 \ / X \ / 1 \ / 1 \ /I 

y' \ = \ -1 \ \ y \ +fi{x + z)\ \ +ez\ \ + Q{x + zf \ 
z' ) \ O) \z ) \-l) \-l) \-l 

where x = 2V + V" , y = V , z = V + V" , and e > and /i are small. 

Next we compute the normal form of (A.6). By finding an appropriate change of vari- 
ables of the form 

(x, y, z) = {Id + fiTi + eT2 + Q){x, y, z) 

where the T are linear operators and Q a quadratic vector field, we obtain after a lengthy 
but straightforward computation 

X =y + fix/2 + Tiz + 0(3) 
(A.7) y' = -x + fiy/2 + Tyz + 0(3) 

z' = -{e + fi)z -T{z^ + (x2 + f)/2) + 0(3). 

where 0{3) denotes terms cubic order or higher terms. Finally, using polar coordinates 
r^ = x^ + y'^ and 6 = cos~^ (^), and dropping all resulting 0(3) terms we obtain the 



normal form 



r = u — \- Lrz 
^2 



(^•') z' = -(e + /.)z-r(z2 + ^) 

9' = -l 

valid in a neighborhood of (e, /i) = (0, 0). 

The dynamics of the normal form equation (A. 8) is quite easy to describe. In the 
corresponding three dimensional phase plane we find two fixed points in the r-z plane with 
r = 0: namely, {r, z) = (0,0) and {r,z) = (0, — (e + ;u)/r). By definition, these correspond 
to the steady solutions U± of the original system (A. 2). The Jacobian matrix J at these 
points are, respectively, 

J(0,0)=f'^/2 / ,landjfo,-(^ + ^)^^-^-^/2 



-{e + fi) J V' r J \ e + fi 

Furthermore, nontrivial stationary solutions of (A. 8) in the r-z plane exist only if fj? + 2fie 
is positive, in which case there is only one possible stationary point P given by 



<-)- V^^ 




A APPENDIX: SURVEY OF EXISTENCE THEORY FOR GENERALIZED KS 56 

which is a non-trivial stationary solution of (A. 8) corresponding to a periodic solution to 
the original system (A. 2). We readily find that the trace and determinant of the Jacobian 
matrix Jp are 

tr(Jp) = -e, det(Jp) = ^^^±^ 

so that the periodic point P undergoes a Hopf bifurcation (in the r-z plane) at e = 0, which 
corresponds precisely to the Kuramoto-Sivashinsky equation. To determine whether it is 
a sub or supercritical bifurcation, one would need to compute higher order terms in the 
normal form, which is beyond the scope of this paper. These periodic orbits emerging from 
P in the r-z plane correspond to quasiperiodic solutions of the original system (A. 2). 

Finally, note that in [CD] a similar normal form was derived in the case when e = 0, 
corresponding to the "classic" KS equation of the form 

r = —fir/2 + 2fz, z = jiz — 2z — 4f , 6' = 1 

from which we can recover (A. 8) in this case by a simple rescaling. In both cases, if 
det (Jp) = ^ 2 ^^ 7^ 0, a straightforward application of the implicit function theorem shows 
that the point P, corresponding to a periodic wave train of (A.l), persists under higher 
order perturbations. 

We further note that in [CD] a full family of periodic solutions to (A. 8) for the Poincare 
return map around the point P and ending with a solitary wave was found. However, it 
was not proved that such a family of quasi periodic solutions persists under higher order 
perturbations. This is not surprising since we are precisely at the Hopf bifurcation point 
and e is the additional parameter needed in KS to carry out a codimension 2 bifurcation 
analysis. In contrast when < \e\ <^ 1, there is a selection of the periodic orbit of the 
Poincare return map and this structure persists for higher order perturbations. 

Let us mention that the full bifurcation analysis for KS is far more complicated: in- 
deed Kent and Elgin [KE] proved the occurrence of a Shi'lnikov bifurcation which leads to 
cascades of period doubling, period multiplying /c-bifurcations and oscillatory homoclinic 
as period is increased. The computation of the bifurcation diagram was also investigated 
numerically for KS [BKJ] and work is still in progress to carry out a similar program for 
gKS by using AUTO continuation software. 

A. 3 The KdV limit: singular perturbation analysis 

We finally conclude our survey of relevant existence results for periodic traveling wave 
solutions of (A.l) by considering a particular singular limit arising in applications to pattern 
formation analysis. In this context, it is usually assumed that e = 1 and < (5 ^ 1 in 
which case (A.l) can be treated as a singular perturbation of the integrable KdV equation 

(A. 9) ut + uux + Uxxx = 0. 
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As is well known, (A. 9) admits a four parameter family of periodic traveling wave solutions 
of the form 

(A. 10) UKdv{x - c{uo,K,p)t;xo,uo,K,p) = uq + 12p^K^cn^ {k{x - xq - c{uo,K,p)t),p) 

with wave speed c(no, K,p) := Sn'^p'^ — 4k^ + uq, where cn(-,p) denotes the standard Jacobi 
elliptic cosine function with elliptic modulus p G [0, 1), and xq, uq, and k are arbitrary real 
numbers corresponding, respectively, to the translation, Galilean, and scaling symmetries 
of the KdV equation (A. 9). 

Here, we are interested in the continuation of these explicit solutions to the KdV equation 
in the singular limit 5 — )■ 0^ in (A.l). In this limit, it was shown in [EMR], by using Fenichel 
theory, that periodic solutions to (A.l) remain close to the above elliptic function solutions 
of the KdV and, in particular, an expansion of these solutions with respect to 5 was obtained. 
Furthermore, we point out that in [BaN] a formal spectral analysis of these perturbed KdV 
waves was conducted. We now briefly describe the associated expansions of the periodic 
traveling waves in the limit 6 — )■ 0"^; see [EMR, BaN, NR2, JNRZ4] for more details. 

Without loss of generality, we can restrict to zero- mean solutions of (A.l). We seek an 
expansion of the associated periodic traveling wave solution of the form 

U\x, t) = UKAvix, t; xo, u*o, K,p) + 6U{x, t) + 0{\6\'^) 

where u* is chosen so that the mean of C/Kdv over one spatial period vanishes. Notice at 
order 5^ the parameters p and k are completely arbitrary, yielding a three parameter family 
of periodic traveling wave solutions of (A. 9) with zero mean. When < (5 <C 1 however, 
we expect only a two parameter family to persist, hence we expect a selection principle to 
manifest itself at the next order between the parameters p and k. 

At the next order one finds that the first-order correction U must satisfy the equation 

(A.ll) H^U'" + {GUKdvU - C{ul, K,p)Uy - CiU{^^y + KU^^y + k^'uZy = 0, 

where ci denotes the first order (5-correction to the wave speed of U . As the linear operator 

2 d^ d 



^■^'^ ^ + ^((Sf^KdV -c(uo,K,p)).) 



has Fredholm index and the kernel of its adjoint operator is spanned by constant functions 
and UKdv{''-,Xf),UQ,K,p)}, it follows that (A.ll) will have a solution precisely when the 
solvability condition 

(A.12) {{uKAyf) = ^\{uKAy^ 

is satisfied. Condition (A.12) yields an explicit selection principle between the elliptic 
modulus p and the scaling parameter k and, as a result, we find for < (5 ^ 1 a two- 
parameter family of periodic traveling wave solutions of (A.l) parameterized by translation 
and the elliptic modulus p or, equivalently, by translation and period. 
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Coming back to the more general case where we ahow UKdv to have non-zero mean, 
we have obtained for < 6 <^ 1 a three dimensional manifold of periodic traveling wave 
solutions of (A.l) parameterized by translation, period, and spatial mean over a period. 

The numerical tests carried out here and in [CDK] suggests that there are finite limit for 
lower and upper bounds of stability as 5 — )■ 0, in agreement with results predicted by formal 
singular perturbation analysis in [BaN] . The purpose of the recent work [ JNRZ4] , although 
still finally relying on elliptic integrals numerical computations of [BaN], is precisely to go 
a step further towards a complete analytic proof of these observations. 

B Appendix: The Swift-Hohenberg equation 

In this appendix, we demonstrate a streamlined proof that spectrally modulationally stable 
periodic traveling wave solutions of the Swift-Hohenberg equation (1.3) are nonlinearly 
stable to small localized perturbations. While this problem has been previously solved 
by Schneider in [Sc] by using a combination of weighted energy estimates, renormalization 
theory, and ingenious nonlinear cancellation technique, all carried out in the Bloch frequency 
domain, our nonlinear analysis rather relies on spatial domain techniques developed in 
[0Z4, JZl, JZ3] in the context of systems of viscous conservation laws. We also carry 
out a numerical spectral stability analysis, demonstrating generality of our techniques. In 
particular, in a specific parameter regime we obtain nice agreement with stability curves 
found by Mielke [Ml]. 

B.l Setup and main result 

The traveling wave solutions of (1.3) are stationary solutions of the PDE 

(B.l) dtu - cd^u + (1 + dlfu -ru + f{u) = 

for some wave speed c G M, i.e. they are solutions of the traveling wave ODE 

(B.2) -cu + u"" + 2u" + (1 - r)u + f{u) = 0. 

Due to the presence of the non-conservative terms, this equation can not be further inte- 
grated, hence the orbits of (B.2) lie in the phase space M^. In particular, it follows that 
periodic orbits u of (B.2) correspond to values {b,c,X) G M®, where X,c ^ M. denote 
the period and speed, respectively, and the vector h = (fci, 62, ^s, ^4) denotes the values of 
(u, u\ u'\ u'") at X = such that 

{u,u ,u",u"'){X;b,c) = b, 

where {u,u',u" ,u"'){-;b,c) is the unique solution of (B.2) so that {u, u' , u" , u"'){0; b, c) = b. 
As usual, we make the following technical assumptions: 

(HI') f eC^, K>5. 
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(H2') The map H : R^ -^ R'^ taking {b,c,X) ^ {u,u' ,u" ,u"'){X;b,c) - 6 is full rank at 
{b,c,X). 

By the Implicit Function Theorem, conditions (H1')-(H2') imply that the set of pe- 
riodic solutions of (1.3) in the vicinity of u, with parameters {b,c,X), forms a smooth 
2-dimensional manifold 

{ (x, t)^U{x-a- c{(3)t; /3) | (q, /3) G R x X } , with X C M. 

In particular, in contrast to the Kuramoto-Sivashinsky equation (1.2), the lack of conser- 
vative structure implies a loss in dimension of the periodic solution manifold. As we will 
see, this has the effect that there are no longer "enough" periodic orbits around to make 
variations in wave speed an admissible perturbation, hence the corresponding linearization 
does not have a (co-periodic) Jordan block at the origin. 

Remark B.l. As noted in [JZN], transversality, (H2'), is necessary for our notion of 
spectral stability hence there is no loss of generality in making this assumption. 

To begin our stability analysis we consider the linearization of (B.l) about a fixed X- 
periodic traveling wave solution u = u{-;b,c). Here we assume without loss of generality 
X = 1. To this end, consider a nearby solution of (B.l) of the form u{x) + v{x,t) with v 
small. Directly substituting this into (B.l) and neglecting quadratic order terms in v leads 
us to the linearized equation {dt — L)v = with 

(B.3) Lv := cvx - v^xxx - 2wxz + {r - l)v - f'{u)v. 

Introducing the one-parameter family of Bloch operators 

L^ := e-^«^Le*«^ ^e ["Vr, ^) 

operating on Lpgj,([0, 1]) parameterizes the spectrum of L as 

0-L2(M)(^) = U O-i2_^,([0,l]) i^d ■ 

ge[-7r,7r) 

We assume the following spectral stability conditions: 

(Dl') a{L) C {A G C|Re(A) < 0} U {0}. 

(D2') (j{L(.) C {A G C|Re(A) < -9\i\'^], for some 6* > and any ^ G [-7r,7r). 

(D3') A = is a simple eigenvalue of Lq. 

By standard spectral perturbation theory [K], (D3') implies that the critical eigenvalue A(.^) 
bifurcating from A = at ^ = is analytic in ^. In particular, under assumption (D2') and 
the symmetry A(^) = A(— ^), the A(^) admits an expansion as ^ — )■ of the form 

(B.4) \{i) = -iai-bf + 0{\i\^), aGM, 6 > 0. 
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Remark B.2. More generally, assume that a real valued, periodic coefficient differential 
operator C whose essential spectrum near the origin admits a Bloch representation as a 
single analytic curve of the form A(^) = "^j^^cxj^^ . Then the symmetry A(^) = A(— .^) 
implies that a2j G M and a2j+i G iM for each j G N. Finally, we note that this observation 
is related to properties of the associated Whitham modulation equations as well; see [JNRZl] 
for details. 

With these preparations in hand, we now state the main theorem of this appendix. 

Theorem B.l. Let u he any steady 1-periodic solution of (B.l) such that (H1')-(H2') and 
(D1')-(D3') hold. Then there exist constants eo > and (7 > such that for any uq with 
11^0 — ^llLi(iR)nH^(R) — ^0; where K is as in assumption (HI'), there exist u a solution of 
(B.l) satisfying u{-,0) = uq and a function '4){-,t) G VF^'°^(M) such that for all t > and 
2 < p < oo we have the estimates 

||n(- + V'(-,t),i)-^(-)llLP(R)<C(l + t)"^('"^)"^||n(-,0)-t2(-)ILi(M)ni^x(R), 
t,^x)(-,t)|lLP(K) <C(l + t) n ''^ MI^(-,0)-u(-)||ii(R)n//X(R), 

3 
M- + ^i-, t),t) - u{-)\\hK(R) <C{l + ty^ \\u{-, 0) - u{-)\\L^(R)nHK{R) , 

3 

Moreover, we have the L^(M.) n H (M.) — )■ L^{M.), 2 < p < oo, nonlinear stability estimate 

(B.6) ||n(-,t)-n(-)||iP(M), IIV'(-, t)llLp{R) < ^ (1 + t)"^('"^) ||n(-, 0) - t2(-)ILi(M)ni/^(R) 

valid for all t > 0. 

The outline of the proof is as follows. Given a fixed periodic solution u{x) of the traveling 
wave ODE (B.2), nearby solutions of the traveling SH equation (B.l) u are investigated by 
defining a nonlinear perturbation variable 

(B.7) v{x, t) := u{x + ip{x, t),t) — u{x), 

where ^ : M x R+ — t- M is a spatial-temporal phase modulation tp to be chosen later such 
that V'(")0) = 0) corresponding to perturbations which are initially spatially localized. By 
a direct calculation, we find that (B.l) is then written 

{dt-L){v-u.,^)=M[v,tl;] 

where L is the linearized operator defined in (B.3) and M denotes a nonlinear remainder 
term depending on derivatives oiv and ip: see (B.15) below. By an application of Duhamel's 
principle, it follows that (B.l) reads 

(B.8) vi-,t) - u^ij{-,t) = e^^vi-,0) + [ e^^^~'^M[v,'ilj]i-,s)ds. 

Jo 
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Of course, we need then a detailed study of the hnearized solution operator e . In the 
present periodic context, this is complicated by the fact that the spectrum of L agrees 
with the essential spectrum and contains a spectral curve touching the imaginary axis. We 
are thus led to separate the exponentially-stable spectrum from this low-Floquet critical 
curve for which, at best, one should only expect polynomial decay of small perturbations. 
A full separation at the linear level would lead to a separation of the nonlinear equation 
into a decay-critical equation and an exponentially slaved equation forced by nonlinear 
coupling terms. Although we perform the linear separation only in an approximate way, 
the purpose of the introduction of ip is to ensure that ipUx contains at the linear level the 
main contribution of the algebraically decaying part of the semigroup, decay of v being 
essentially slaved. 

B.2 Linearized estimates 

We begin our analysis deriving decay rates for the semigroup e of the linearized operator 
L defined in (B.3). With the inverse Bloch transform representation (1.13) of the solution 
operator in mind, we begin by making more precise the statements about the critical curve. 
Since is separated from the rest of the spectrum of Lq, we first note that by standard 
spectral perturbation theory [K], assumption (D3') implies that there exist ^o G]0, 7r[, eo > 0, 
an analytic curve A : [— ■^o^'^o] — ^ ^(0,£o) such that, when |^| < ,^0) ^(-^5) l^ ^(0,eo) = 
{A(^)}. and, when |^| < ^o, dual right and left eigenfunctions q{(,,-) and q{C,-) of Lg 
associated with A(.^), analytic in ^ and such that q{0,-) = Ux- We may thus define the 
critical spectral projection -P(^) through 

Symmetry of the spectrum and (D2') imply that A satisfies (B.4). 

Now, we choose a smooth cutoff function (p : [— tt, vr) — )■ [0, 1] such that 

, ... _ r 1, if 1^1 < ^0/2 

and decompose the linearized solution operator S{t) = e^* as S{t) = S^{t) + S^^{t), with 
S'it)gix) := r e*«"0(O[P(e)e^«*5(C,-)](x)rfe 



e 



^^^^sr/^V^CS)* 



0(e)e"^^^*g(e, x)m, ■),m, ■))Li^Am)d^ 



S''{t)g{x) := / e^«-[(l - 0(0^(0) e'^^'dit ■)]{x)dC. 



— TV 



Note, by the way, that the Green kernel 

G'ix,t;y):=S'it)5yix) 
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is given by 



G'{x,t;y)= I e^«"-^)0(e)e^(«)*g(e,x)g(^,y)*d^; 



see the proof of Lemma 3.2 for some details. General considerations about semigroups lead 
for S to the following bounds, whose proof is essentially the same as Proposition 3.1. 

Proposition B.2. Under assumptions (H1')-(H2') and (D1')-(D3'), there exist constants 
C,6 > such that for all 2 < p < oo, < Ah + k = k < K + I, < Ami + m2+m^<K, 
K as in (HI '), s = 0, 1 and t > 



di'di^s''{t)di^g 



\d^id^2sii(t)dT^g\ 



H" 



LP( 



<Ct 4 e ''^Mhs 



1/1 l'\ 4mi+m2+m3 



< Ct 4V2 P 



W^ (*)5'ILp(r) - ^^ *ll5ll//i(iR)- 

Aiming at identifying up to remainders terms of form S {t)g with terms of form Uxh, 
we decompose further by expanding q{^,-) = Ux + 0(0 ^^^ thus write G^{x,t;y) = 
Ux{x)e{x, t; y) + & {x, t; y) with 



Since 



e(x,t;y) 
&{x,t-y) 

e{x,t;y)g{y)dy = 



e^^^--y^{i)e^^^)'q{Urdi 

TV 

' e*?(--?/)0(^)e^«)* (^(^, a;) - qiO, x)) qi^, yYd^ 

TV 



&{x, t; y)dy = / e^«-0(e)e^(«)* {q{i, x) - (?(0, x)) m, ^.^ •))L2,.([0,i])dC, 

a direct application of Hausdorff- Young inequality (1.15) yields the following bounds. 

Proposition B.3. Under the assumptions (H1')-(H2') and (D1)-(D3), the low-frequency 
Green function G^{x,t;y) of (B.3) may be decomposed as 

G^{x, t; y) = u{x)e{x, t; y) + G^{x, t; y) 

where for allt>0 and l<q<2<p< oo, < j,l,s,j + 1 < K + 1, the residual G^ {x,t; y) 
satisfies 

(B.9) 



didldlG'{;t;y)g{y)dy 



_ 1 / i_i\ i_i 

<C(l+t) ^9 v) 2 2 



mm 



LP( 



Furthermore, for all t >0, l<g<2<p< oo, < j,l,s,j + I < K + l,we have 



(B.IO) 



didtdyei-,t;y) g{y)dy 



1 / 1 1 A (j+0 

<G{l+t) H9 p) 2 



mi'ii 



LP( 
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Finally, we combine the various exponential and algebraic bounds derived above to 
obtain decay estimates on the Green function 



(B.ll) 



Gix,t;y) = Sit)6yix) 



associated with the full solution operator S{t) = e . To this end, let x be a smooth real 
valued cutoff function defined on [0, oo) such that < x ^ 1 and 



x{t) 



0, if < t < 1 

1, if t > 2 



and define 



e{x,t;y) := x{t)e{x,t;y). 
Using the Hausdorff- Young inequality (1-15) and the triangle inequality 



F{-,t,y)g{y)dy 



< C'll5'llLi(K)Sup||i^(-,i,; 



iLPfR) ' 



LP( 



we obtain the following estimates. 

Corollary B.4. Under assumptions (H1')-(H2') and (D1')-(D3'), the Green function 
G{x,t;y) of (B.ll) decomposes as 

G{x, t; y) = u{x)e{x, t; y) + G{x, t; y) 

where, for some C,9 > and all t > 0, 2 < p < oo, and 0<j,l,s,j + l<K + l we have 



(B.12) 

(B.13) 
(B.14) 



d'G{;t;y)g{y)dy 



G{-,t;y)g{y)dy 



didld'e{-,t;y)g{y)dy 



1/1 1 \ s 



<Ct ^V2 pJ ■*(! + *) H2 l) 2 + 4 



IS'IIli 



LP( 



LP( 



LP( 



<C(l+t) 



-i 1-i - 



pJ 2 



l5llLi{]R)n//i 



if, i\ (j+fc) 
<C(l + t) H' p) 2 



\9\\m 



and e{x,t;y) = for < t < 1 . 

B.3 Nonlinear analysis 

Given the linearized bounds on the linearized solution operator S{t) = e^* derived in the 
previous section, we are now in position to consider the effect of the small nonlinear terms 
that were omitted in obtaining the linearized equation. For this purpose, we introduce 
v{x, t) = u{x + 4>{x, t),t) — u{x) as in (B.7). Direct calculations similar to those detailed in 
the proof of Lemma 3.3 provide the following lemma. 
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Lemma B.5. For v as above, the equation is 

(B.15) idt-L){v-'il^u^)=M, M=Q + d^n + dtS + r 

where 

Q = - {f{u + v)- f{u) - df{u)v) 

n = i)tv- 2—-—^v^ - ——^dcc ( , , , ^x ( -:——-Vx ]] - dxi r—r^d, 



1 + ^^"" 1 +^2;"" VI +^2:"" Vl + ^^""yy "" Vl + V'a;''^ Vl +V'a; ^ 



~Uxxx Cx I :; ; ; dx I z ; ; Ux I I Ox I :; ; ; Uxx I C™ 



-, , , '-I'XXX ^X \ -, , , ^X 1-1, , '-I'X I I ^x \ -, , , f^xx I '-'x I -I , / "'X 

i + ipx Vi + V'x Vi + V'x J J \i + i>x J ""xi + Tpx 

S = -vipx 

r = -V'x ((1 - r)v + /(n + v)- f{v)) . 

By Duhamel's principle, we obtain the announced implicit representation (B.8), so that 
we may express the phase ^ and the nonlinear perturbation variable v implicitly as 

(B.16) ^{x,t) = - I e{x,t;y)v{y,0)dy- / / e{x,t - s;y)Miy,s)dyds 

JR Jo Jr 

and 

(B.17) v{x, t)= [ G{x, t; y)v{y, 0)dy + f f G{x, t - s; y)Af{y, s)dy ds. 

JR Jo Jr 

Furthermore, recalling that e{x,s;y) = for < s < 1, we find by difi'erentiating (B.16) 
that 

(B.18) d{d^4^{x, t) = - f d{d^e{x, t; y)v{y, 0)dy -If a/a^e(x, t - s; y)M{y, s)dy ds. 

Jm. Jo Jr 

for < J < 1 and < k < K + 1. To apply a standard contraction-mapping argument and 
solve locally in time our Cauchy problem (B.17)-(B.18) with {v, ipt, ipx) G H^{R) x H^(R) x 
H^^^(K), we only need besides our linear bounds the following nonlinear damping energy 
estimate whose proof is entirely similar to the one of Proposition 3.4. 

Proposition B.6. Assuming (HI '), there exist positive constants 9, C and Eq such that if 
V and tp solve (B.15) on [0,T] for some T > and 

sup \\{v,iJx){t)\\HKm+ sup \\iJt{t)\\HK-i(^) < eo 

te[0,T] t£[0,T] 

then, for all <t <T, 



Ht)fHK(^.<Ce-''\\vml^, 



(B.19) /•* . X / \ 

+ C / e-^(*-^) (ll^(^)lli2(K) + ||V.(s)ll Whm) + \\Ms)\\Ik-.^^)) ds. 
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With the above preparations in hand, we are now prepared to state the main technical 
lemma leading to the proof of Theorem B.l. To this end, associated with the solution 
{u,^t,'<l^x) of the integral system (B.17) and (B.18) considered in the previous section, 
define 

(B.20) 7i{t):= sup \\{v,iJtAx,^xx){s)\\HKo^Al + sf'\ 

0<s<t 

By standard short-time /7^(M) existence theory, the function t] is continuous so long as it 
remains sufficiently small. Using the linearized estimates of Section B.2 we now prove that 
if r]{0) is sufficiently small then ij^t) remains small for all t > 0. 

Lemma B.3. Under assumptions (H1')-(H2') and (D1')-(D3'), there exist positive con- 
stants C and £0 such that if v{0) is such that for some T > 

Eo ■■= \\viO)\\L^R)nHK(K) < £o and r]{T) < eq 

then, for all < t < T , 

(B.21) 7]{t)<C{Eo + r]{tf). 

Proof. First, note that under the above smallness assumption 

\mmLHR)nLm<Cv{t?{l + t)-'/'. 

Applying now the bounds of Corollary B.4 to representations (B.17)-(B.18), we obtain for 
any 2 < p < oo the bounds 

lk(t)|lL.(M)<C(l + t)'Ki-^)-^i5;o 

(B.22) + C7]{tf f {t- s)~K^'p) (1 + t - s)-Ki~p)~^ (1 + syl^s 

Jo 



<C(Eo + vitf){l + t) H^ p) ^ 



and 



1/1 i\ 1 



ll(V'i,V'x)||vKi^+i,P(M)<C(l + *) '^' ^^ 'Eo 



(B.23) + C7]{tf I (1 + t - s)-5(i-i/p)-^ (1 + s)--^/2(is 

Jo 

< c {Eo + 7]{tf) {1 + ty'^^^'^p)'K 

Finally, since the smallness assumption guarantees that we can apply Proposition (B.6) it 
follows that 

ft 



\V(t}\\]jK( 



<Ce-''Ei + C{Eo + 7^{tff f e-'('-'\l + s)-^/^ds 

Jo 

< Ce-^'El + C{E^ + r^itfY (1 + ty^^ 



\-3/2 



which, together with (B.23) , completes the proof. D 
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Assuming without loss of generality that the constant C in (B.21) is larger than 1, we 
obtain, since rj is continuous and 7/(0) = ||t'(0)||//xm) < Eq, by continuous induction that if 
AC'^Eo < 1 then r?(t) < 2CEo for all t > 0. Recahing (B.22), (B.23) and the definition of r/ 
yields (B.5). 

Using Corollary B.4 again gives for 2 < p < oo 

LP(R) <CEo{l + ty'^^^'^) +CEl j (l + t-s)"K^-p)(i + s)-ids 



<CSo(l + t) ^^^ ^ • 
Now by definition we have that 

u{x, t) — u{x) = v{x, t) + {u{x, t) — u{x + ip{x, t),t)) , 

hence 

_ifi_i 

\\u{t)-u\\LP(R)<\\v{t)\\Lp(^) + \\Ux\\L'^([0,l])\\'>P{t)\\LP{R)<CEo{l + t) H P . 

This completes the proof of Theorem B.l, establishing the nonlinear L^ n H — )■ L^ asymp- 
totic stability of the underlying periodic traveling wave u under the structural and spectral 
assumptions (H1')-(H2') and (D1')-(D3'). 

B.4 Numerical stability analysis 

We demonstrate the numerical stability verification method described in Section 2.1.3 in 
the context of the Swift-Hohenberg equation. To match with [Sc], we choose f{u) = u^, set 
£ = \fr and rewrite (1.3) as 

(B.24) dtu = (e^ - \)u - 2dlu - d*u - u^. 

Stationary traveling wave solutions of (B.24) satisfy, 

(B.25) (e^ - l)u - 2u" - u"" - u^ = 0. 

When £ is small there is a three parameter family of stationary periodic solutions approxi- 
mated by 



(B.26) uo(a;, 0, e)\x\ = 2Re(e(\/l - 4a;2/\/3)e^( 



l+eu})x i4>\ 



and with period X = 27r/(l + ew); see [Sc]. About n, a periodic stationary wave solution 
of (B.24), the linearized evolution is described by [dt — L) = where 

(B.27) Lu := (e — l)u — 2uxx — Uxxxx — Sti u, 

Note that the fact we deal with stationary solutions make L and Lq self-adjoint. In 
particular, assumption (Dl') of the previous section is a consequence of (D2') and the 



B APPENDIX: THE SWIFT-HOHENBERG EQUATION 67 

spectrum of L lies in M so that any transition to stability/instability must be marked 
by an eigenvalue passing through the origin for some possibly non-zero Bloch frequency. 
Furthermore, in [Eck, Ml, MS] it was analytically verified that for e > small in (B.24) 
the solutions approximated by uo{oj,4',e) are spectrally stable, i.e. (Tj^2 (ig) (L) C (— oo,0], 
provided that 

(B.28) \4oj^\<\ + 0{e). 

Below, we apply an appropriate modification of the numerical protocol introduced in Section 
2.1.3 above to analyze the stability of such small amplitude periodic wave trains of (B.24). 
Modifying adequately the proof of Lemma 2.2, we find that (Tj^2 (jj) (L) C (— oo,e^]. We 
first check that there is only one spectral branch of the associated Bloch operators L^ 
bifurcating from the origin at ^ = 0, and then track its location. In this case, the Evans 
function D{X, S,) can be expanded to second order as 

D{X, C) = aoX + ai^ + asA^ + asX^ + a^f + O {\X\^ + |e|^) . 

Assuming that indeed only one critical spectral branch bifurcates from the (A, ^) = (0, 0) 
state, expanding the critical spectral branch as A(^) = a^ + /3^^ + 0(|^p) it follows that 
the coefficients a and /3 are related to the above Evans function expansion via the formulas 

oi , ^ a2a^ + 030 + 04 

a = and p = . 

oo «o 

Recalling Remark B.2 together with the fact that here the spectrum is real, we note that 
a = so that /? = — — , and thus there are only two coefficients of the Taylor expansion 
of the Evans function that need to be computed in order to check stability for |(A, ^)| <^ 1 
with stability in this small-Floquet regime corresponding to the condition that /3 < 0. 

For e = 0.187, uj = 0.1, and X = 6.1678 we compute the Evans function on a contour 
d{B{0,R)/B{0,ro) n {A|Re(A) > 0}) where R = e^ and tq = IQ-^ using 1001 Floquet 
parameters, ^, evenly spaced in the interval [—tt/X,tt/X], and an adaptive mesh in the A 
contour requiring relative error be less than 0.2 between consecutive contour points. We 
then use the Taylor coefficients method to find /3 = — 3.874-|-0.000i via a convergence study 
requiring convergence in relative error between successive iterations of /3 be less than 0.01. 
As a check on the accuracy of our computation, we note that a = 5.2664e — 11 + 1.9970e — 
07i ~ 0. Next we determine the maximum modulus root of the Evans function for |^| = 0.2, 
^ € C using a A contour of radius 2, yielding the bound, max(A(^)) = 1.1099. 

We find that fco = 1.7420e — 4 suffices for breaking the Floquet parameter interval into 
small and large modulus values. For |^| > /cq we compute the Evans function on a semicircle 
of radius 0.1 passing through the origin requiring that the relative error between consecutive 
contour points not exceed 0.2. Then we compute the Evans function on the same semicircle 
shifted left by le — 4, this time taking |^| < /cq. In all computations, we find the winding 
number is consistent with spectral stability. 
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We also find the stability boundary, though this time we do not use Taylor coefficients 
because the algorithm for a and /3 becomes numerically poor conditioned since a becomes 
small near much of the stability boundary; see Remark 2.4 above. To circumvent this 
difficulty, we find the stability boundary as given by computing the Evans function on a 
contour 5(5(0, e^)/i?(0, ro)n {A I Re(A) > 0}), ro < e^, to determine the presence of spectra 
as uj varies, and then take the limit of the location of the stability boundary found in this 
way as ro — )■ 0. We take care to choose via a convergence study a sufficiently tight tolerance 
setting in the RKF integration routine to obtain accurate results. We plot the stability 
boundary and a comparison with the analytical bound given by Eckhaus, Schneider, and 
Mielke [Eck, Ml, Sc] in Figure 10. 

Finally, we verify numerically the eigenvalue picture given in Figure 1 of [Sc] by com- 
puting the Evans function on a circle and then solving for the roots using the method of 
moments as described in Section 2.1.5; see Figure 11. 
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Figure 10: In (a) we plot the stability boundary for the Swift-Hohenberg equation, and in 
(b) we compare the numerical stability boundary with the analytic one in (B.28), plotted 
here using dashed lines, found by Eckhaus, Schneider, and Mielke [Ml, MS]. This provides 
a check on our numerics against a well known result. The analytic curve is found by 
rigorous perturbation analysis in [Eck, Ml, MS]. Our numerics give a nice extension to 
large amplitudes in the non-perturbative regime. We find the stability boundary with 
relative error tolerance of 0.01. In (a) we do a convergence study as ro — )■ as described 
previously. In (b) we compute the Evans function on the contour 0(5(0, e^)/i?(0, ro) n 
{A I Re(A) > 0}), ro = le — 8, with tolerance set at le-12 in the RKF routine. 



C Appendix: Computational statistics 

In this short appendix we detail computational statistics for some typical values in the 
numerical studies reported in previous sections. All numerical computations were performed 
on either a Mac Pro with 2 Quad-Core Intel Xeon processors with speed 2.26 GHz, the super 
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Figure 11: Plot of two smallest eigenvalues of the Swift-Hotienberg system against the 
Floquet parameter ^ for e = 0.05, oj = 0.28858; cf. Figure 1 in [Sc]. 



computer Quarry^® at Indiana University, an Intel(R) Pentium(R) 4 CPU with speed 2.80 
GHz running Windows XP, or a MacBook with an Intel Core 2 Duo processor with speed 
2.0 GHz. 

Carrying out the stability study for the classical Kuramoto-Sivashinsky equation with 
q = 6 and X = 6.3, corresponding to the q = 6 row of Table 1, with the MacPro using 8- 
cores via MATLAB's parallel computing toolbox, it took 11.46 seconds to perform the High 
Frequency Evans function computation, 159.0 seconds to compute the Taylor coefficients 
{aj}j=i^2 and {I3j}j=i^2, 26.85 seconds to compute max^ |Aj(^)|, 169.5 seconds to compute 
the Evans function when^^ fco < |.^| < vr/X, and 40.89 sec to compute the Evans function 
for 1^1 < /cq. For these same parameters computing SpectrUW on the Windows machine 
using the Maple kernel took 26 seconds using 5 Fourier modes and 30 Floquet parameters. 
It took 55 seconds using 10 Fourier modes and 30 Floquet parameters. 

The stability study computation times for the Swifth-Hohenberg equation, those dis- 
cussed in Appendix B.4 using the MacPro with 8-cores are as follows. The High Frequency 
Evans function computation took 5.638 seconds, the Taylor coefficients took 81.73 sec- 
onds to compute, finding max |Aj(i^)| took 1.110 seconds, computing the Evans function for 
ko < \^\ < tt/X took 11.30 seconds, and computing the Evans function for |^| < /cq took 
8.589 seconds. 

The stability diagram in Figure 6 for the full parameter study of the generalized Kuramoto- 
Sivashinsky equation took roughly 6 days on 6 nodes with 8 cores each on Quarry using 
MATLAB'S parallel computing capabilities. We took steps sizes of 0.25 in e and 0.1 in 
period X to check for stability. Furthermore, the stability picture in Figure 10 for the 
Swift-Hohenberg equation took roughly 3 days running on the MacPro. Finally, concerning 
the time evolution studies in Section 2.3, those over time intervals of less than 50 took 
between 5-30 minutes each to run while those on larger time intervals took approximately 



^^This is a cluster composed of IBM HS21 Bladeservers and IBM iDataPlex dx340 rack-mounted servers. 
^'^Here, we are following Section 2.1.3. In particular, kg is the low-frequency cutoff in Step 2(b) of the 



numerical procedure outlined in Section 2.1.3. 
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5 hours each to complete on the Macbook described at the beginning of this section. 

D Appendix: Numerical Evans algorithm 

Here, we briefly describe the numerical algorithm on which our Evans function computations 
are based, and explain the terminology (scaled lifted polar coordinate method; unsealed 
lifted polar coordinate method) used elsewhere in the paper. Recall the standard definition 
(2.13) of the periodic Evans function as D{^, A) := det($(X) — e*^"'''!), where ^(x. A) is the 
monodromy matrix of eigenvalue ODE (2.12), that is, a matrix- valued solution with initial 
condition ^(0, A) = I. 

As pointed out in [0Z2], this is deceptively simple to code as compared to the homoclinic 
or front-type Evans function defined on the whole line [AGJ, PW, GZ], in the sense that 
there are no issues with infinite domains and imposition of asymptotic boundary conditions 
at X = iboo, so that one may in a matter of moments develop a naive algorithm that is 
serviceable for moderate values of A and X, using any standard ODE evolution algorithm. 

However, as described in [BJNRZl], for the type of global parameter exploration carried 
out here, inherently involving large A and X values, such a naive algorithm is essentially 
useless. This is most easily seen by the asymptotic relation [G2, SS, Z4] between the periodic 
and homoclinic Evans functions as period X — )■ oo, as we now describe. 

D.l Balanced and rescaled Evans functions 

Start with the eigenvalue equation Y'(x; A) = M{x, A)Y(x; A), x G [0, X], Away from a finite 
set of curves Cj{k) determined by the dispersion relation ik G a{M{0,X)), A; G M, where a 
denotes spectrum, the eigenvalues of E[(0, A) have non-vanishing real part. For each j G N, 
set 

Cj := {A G C I A = Cj{k) for some /c G M} 

and let C denotes the compliment in C of the set DjCj. On C, we can define the real- 
valued function n such that n(A) equals the number of eigenvalues of 11(0, A) with negative 
real part noting, in particular, that n is locally constant on each connected component of 
C. Similarly, we define the real valued functions a± such that a-^-{X), respectively a„(A), 
equals the sum of the positive, respectively negative, real part eigenvalues of 11(0, A): as 
above, the functions a± are locally constant on each connected component of C Following 
[G2], we note that, by Abel's formula, D{X,(^) may be written alternatively as 

(D.l) D{\, e) = D{X, C)e~ ^^'^ tr(H(j/,A))dy^ 

where 

(D.2) Z)(A, i) := det(^(0)^(-X/2)-i - e'^^'^{0)^{X/2)-^) 

is a "balanced" periodic Evans function defined symmetrically about x = 0. 
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Definition D.l. For X £ C and S, £ [—Tr/X,iT/X), we define the rescaled balanced periodic 
Evans function as 

(D.3) ^(A, e) := e("- W-"+ WW2e-"W^«^Z)(A, 0- 

The rescaling in (D.3) is designed to cancel the exponential growth in (D.2) with respect 
to the period X of the pieces ^(0)*(-X/2)~^ and '^{0)^{X/2y^ in the constant-coefficient 
case IHI = constant. Indeed, suppose that v,-^ are a sequence of periodic waves with period 
X, converging as X — t- oo to a solitary wave, or homoclinic, solution u°°, and index the as- 
sociated sequence of periodic Evans functions by D^{X, ^). Then, under mild assumptions, 
we have [Z4]: 

(D.4) lim 6^(A,e) = I)hom(A), 

X->oo 

where -Dhom(A) denotes the associated homoclinic Evans function defined in [PW, GZ]. 

Comparing with (D.l) and (D.3), we see that if D-^ is uniformly bounded then D-^ 
and L)-^ exhibit exponential growth with respect to X. Thus, by (D.4), together with 
the well-known fact (see, e.g., [GZ]) that the homoclinic Evans function Dhoiru hence also 
D^ , is uniformly bounded on compact A-domains, that both the usual and balanced Evans 
functions D and D exhibit exponential growth with respect to X as X — )■ oo, hence 
become numerically impractical for large X or for moderate periods X and large A (leading 
to large a±) in the sense that even small variations in A will lead to excessive winding in 
the image I?(A,^) that requires fine resolution to track for purposes of winding number 
computations. 

For our high-frequency winding number estimates, therefore, it is crucial to use the 
rescaled Evans function D in place of D or D. Near A = on the other hand there are 
several curves on which the spectra of ]HI(0) becomes imaginary, across which n, a_, and 
a-i_ undergo discontinuities. This does not occur on {A | ReA > 0} \ {0}, so does not come 
into play in high-frequency estimates; however, for our tracking of spectral curves CjH) near 
zero using analyticity/the method of moments, it is evidently a problem. Thus, for our low- 
frequency computations, we use the unrescaled Evans function, to avoid loss of analyticity 
that would otherwise occur near A = 0, a subtle but important point. 

D.2 The lifted Evans function 

Up to this point we have only discussed winding with respect to A, ignoring questions 
of numerical conditioning. However, these are quite relevant also for large X, given the 
exponential growth already described. Indeed, one faces all of the issues described for 
homoclinic or front-type Evans functions in [Zl], arising from variations in growth rate 
of various scalar modes in the solution operator ^(A,x), in which faster-growing modes 
dominate slower-growing modes through cancellation/loss of significant digits. 

Similarly as in the homoclinic case [Br, BrZ, HuZ], these can be avoided by the device 
introduced in [BJNRZl] of working with "lifted equations" for which the periodic Evans 
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function appears as a Wronskian of certain bases of solutions. Namely, we may recall 
Gardner [Gl], and his simple but important observation 

(D.5) D{^, A) := det(^(X) - e^«^Id) = det (^^[^^ "'^j^^^) 

expressing D as an exterior product of evolving solutions for equations 

Id\ , „ , fe'^^ld\ 



with data ( -. , at a; = and ( , , ] at x = X. 

\idj V Id ; 

Working with (D.5) allows us to apply any of the efficient algorithms that have been de- 
veloped in the homoclinic Evans function setting for evaluation of Wronskians, in particular 
the exterior product method of [Br, BrZ] in which minors of columns in (D.5) are computed 
as evolving exterior products, or the polar coordinate method of [HuZ], in which they are 
computed instead as evolving orthonormal subspace/scalar "radius" pairs; see [HuZ] for 
further details. Both of these methods are by now standard and have been implemented 
in STABLAB, a Matlab-based platform for numerical Evans function computations; see 
[BHZ2]. These methods can thus readily be implemented once we have made the refor- 
mulation (D.5), and the resulting algorithm is numerically well-conditioned by the same 
analysis sed in [Zl] to study the homoclinic case. 

The results of this apparently simple change are dramatic, extending the range of \X\ 
values that can be computed by up to three orders of magnitude; see [BJNRZl]. Indeed, the 
efficiency of the lifted balanced periodic method appears to be quite comparable with that 
of the well-established algorithms used to evaluate the homoclinic Evans function, making 
the numerical Evans function approach a practical alternative to the Galerkin methods that 
have been used in past literature [EST, CKTR, CDK]. 

D.3 Optimization across ^ 

Finally, we mention for the reader who may wish to perform such computations indepen- 
dently a final detail that greatly speeds up computations with varying Bloch frequency ^. 
Namely, rather than computing a new solution of (D.6) for each value of ^, we may, by the 
decoupled nature of the fiow in first and second coordinates, compute the result for a single 

solution initialized as I ^ , ) at a; = X, then multiply the first entry of the resulting solution 



Jd^ 

by e*^ to obtain the desired suite of solutions for varying ^. This means that we need 
in practice only compute solutions of the eigenvalue ODE once for each value of A, with 
variation in ^ introduced by computationally negligible linear algebraic manipulations, for 
a considerable savings. 

See [BJNRZl] for further discussion/comparison of results. 
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E Appendix: Behavior near stability boundaries 

Finally, we consider the difference in behavior expected as we cross a stability boundary 
across which hyperbolicity of the first-order Whitham system is lost vs. behavior expected 
as we cross a stability boundary where strict hyperbolicity is maintained but the second- 
order diffusion coefficient changes sign. 

Introduce a bifurcation parameter ij, with the stability boundary assumed to correspond 
to rj = 0. In the second case, by standard spectral perturbation theory/separation to first 
order of modes, we have Aj(^, ry) = iaj{rj)$, — bj{r])^'^ +icj{'q)£,^ — dj{rj)^'^ + C'(C^), where Oj, 
bj , Cj , dj (by complex symmetry of eigenvalues of real- valued operators) are real and 

bj = -b*r] + o{ri), dj = d^ + o{l), 

so that ReXj{£,,7]) r^ r/6*,^^ — d*^^ and (after a brief calculation) 

(E.l) maxReAj(^, 77) ~ r/ . 

In the first case, we note, rather, that, by the connection to the second-order Whitham 
approximation, Xj{^, rj) = ity^a*^ -|- 6*^^ -|- o(^^), whence 

(E.2) maxReAj(.^, ry) ~ r/. 

We may conclude that the transition involving loss of hyperbolicity is more drastic, featuring 
exponential growth at rate roughly the square root of the rate expected for a comparable 
transition across a boundary involving loss of diffusivity. 
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